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AN ELECTRON-POSITRON JET TODEL FOR THE GALACTIC CENTER 


M. L. Dumis 

NASA/CODDARD SPACE FLIGHT CENTER 
GREEraELT, m 20771 
and 

DEPARTMENT OF PHYSICS and ASTRONOMY 
UNIVERSITY OP MARYLAIH) 
COLLEGE PARK, MD 20742 


ABSTRACT 


High energy observatlono of the galactic center on the 
subparsec scale oeea to be consistent uith electron-positron 
production in the fom of relativistic Jets. These jets could be 
produced by a '-10® black hole dynano transporting pairs euoy from 
the tnasoivc core* An elcctroziagnetlc cascade shower would develop 
first fron aabient soft photons and then non-linearly; the shower 
using itself as a scattering nediun. This is suited to producing, 
cooling and transporting pairs to the observed annihilation region* 
It is possible the center of our galaxy is a nlniature version of 
more powerful active galactic nuclei that e^diibit jet activity* 

OBSERVATIONS 

Interesting activity from the heart of our galaxy has been 
observed over the pr^st few years* Specif ically, a luminous line at 
0*311 MeV has been detected which is ouprtolngly norrow and 
variable*^ In addition, the region from which this line emanates 
coincides with a small and less luminous radio source. ^ A short 
list of these observations follows: 

1. A 10^ object at the galactic center could be inferred fron 
the widths of Neon lines In HII clouds orbiting the central 
region* ^ 

2. Luminosities are: 10^® erg/s in continuum gamma rays**; 10^^ 
erg/s in 0.511 MeV llne^; 10^^ erg/a in hard X-rays**; 10^5 erg/s 
in soft X-rays**; 10**® erg/s in UV^; 10**^ erg/s in IR^ and 10^^ 
erg/s in radio. ^ 

3. The 0.511 KeV line is variable on approximately a six month time 
scale • ^ 

4. The line was measured at (510.9 t 0.25)kcV.^ 

5. A compact radio source is resolved to less than 10^^ cm in 
slze.^ 
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Two laportant conatralnto coae fron the lino* Flrot, the clx 
Qonth variability conatraino the oise of the rcsion producins the 
line to be less than or equal to six light oonths in extent* It is 
poQoible the radio source is an even core otrlngent conotralnt on 
sire. Secondly, the annihilation region cuot be separated froa any 
massive object by at least '^ 10 ^[ 2 Qi/c^J('^ 10 ^S fo^ a 10® black 
hole) so that the line will not be gravitationally redohlftcd cut of 
the observed narrow energy range* Since the compact radio object 
exhibits spectral properties oiniliar to the cores of active 
galaxies, this coupled with the observed higher energy phenomena 
leads one to consider scaled dovm versions of power sources for 
active galaxies* 


MODEL 

A black hole dynamo producing relativistic e*^-e jets of thc*^ 
type proposed by Lovelace',® for double radio sources will be 
considered here* Other beam models have been considered by Brown®, 
Blandford^® and Novikov.!^ Relativistic electron-positron jets seem 
necessary because: 

1. A line is observed at 0*511 MeV implying annihilation of e*^-e" 
pairs* 

2* A jet would beam pairs into an annihilation region far enough 
from the Dassivc engine to be consistent with the width of the 
line and bypascco the problem of producing pairs to close to the 
engine. 

3. The natural propagation of a relativistic jet would convert 
kinetic energy into rest mass of pairs thereby simultaneously 
producing them and transporting them out* 

4. Relativistic beaming maintains a confined interaction region in 
which high energy particles can Interact with one another* 

A dynamo, possibly powering a quasar, produces a luminosity, 

L- 10‘***(M/10®M^)2(b/ 103G)2 erg/s.® For the galactic center a 
^10® Hq black hole and line luminosity of *^10®^ erg/s results in a 
dynamo magnetic field of 30 Gauss* Using this field, an 
accelerating potential can be calculated. It is V--10^®(M/10 ®Hq) 
(B/10®G) V V. The resulting electron flow from the engine la 

I^1037(m/ 10®J1^,)(B/10®G) e/s e/s. These colllnated fast 

electrons will then interact with the dominant scattering 
background. This background Initially appears to be the Infrared 
radiation however the majority of interstellar dust that could be 
correlated with IR is believed to be outside of the central parsec® 
thereby making it impossible for electrons to cool against it* 
However erg/s in ionizing UV must be produced from a 3*2x10**^ 

region within the central parsec Co be consistent with Hll 
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Figure K Scheraatlc of Electron-Positron Jet at Galactic Center 


observations.^ An accretion disc eralttlng W at a density of 
cnT^ from *^10 cm^ could supply a cooling rn'dlum for fast 
electrons* this setting a cross-sectional radius of cm for 

the Jet (Figure 1). Since the compact radio source Is smaller than 
^lO'^ cm* it would seem the hulk of fast particles are not 
propagated past that distance. The radio therefore sets a length of 
the jet at t-lO^^ cm. Knowledge of the observed O.Sll MeV line flux 
and necessary photon reaction rate quite Independent ly enables one 
to constraint the dimensions of the source.'* For beam geometry it 
is found the length and cross-sectional radius irxist be related 
by a*/f-l0** cm for photon-photon Interactions to occur. The UV-ravUo 
estimates of the Jet dimensions are then consistent with constraints 
from phot on -phot on reactions. 



i 


ORIGINAL PAGE (3 
OF POOR QUALITY 


CASCADE SKOUBR 

In the vlnlclty of the dynooo ("*10*2 co) on olectroaagnotlc 
caecado will develop by tho Interaction of colllcatod foot electrone 
with UV. Pairs will be produced, In a linear fashion, aa kinetic 
energy la degraded fron oV to -10* * eV, the pair production 

threshold. At a dlatance of -10*^ ca, pair production froa the UV 
will cease bocauoe the UV density la dccroaalng with Increaelng 
dlatance. The yield of palre , f^, con be calculated from coupled 
relatlvletlc Boltsnann tranoport equations. * ^ The spectrum behaves 
roughly aa a power law. 


Linear: f^(t,E) - (D 

Here c Is the dlatance Into the shower from the Injection point at 
-10*2 cm, Eq la the primary energy from the engine and E la the 
energy. The optical depth of the shower particles to UV Is 

T(C) - 10 fl-(^)’l (2) 

where aQ>lO*2 cm la the engine radius. The flux of shower partlclos 
has then Increased from -10^** e/s Injected by the engine to 
-1038 eVs fia a reoult of linear cascade action* 

Alchou^h pair production froa the UV has ceased, continued 
cooling of pairs from the UV (Figure I, Region I) will produce hard 
photons energetic enough to sustain pair production in a non-linear 
fashion* The non-llnear shower proceeds as follows* Hard photons 
produced by inverse Compton of the UV will pair produce. Those pairs 
then cool by Inverse Compton, producing more hard photons which 
further pair produce and so on. However, since the UV density Is 
dropping with distance and the shower population is growing, at some 
point the soft photons in the shower will dominate the UV (Figure I, 
Region II). The cascade then proceeds using soft photons In the 
shower as a medium until pair production threshold is reached. The 
end-products are then cooled pairs and hard X-rays, below threshold, 
of comparable luminosities (^10'^^ erg/s). The yield of pairs, f^, 
can he calculated for tho non-linear shower from coupled transport 
equations.'** The simpler solutions yield power laws of the form 

Non-Une«r: f^(s.E) - fl - X (^ - 1)1"^ (3) 

- o 

Here z is the distance into the shower, Cq« 10'^ cm is the starting 
position of the non-llnear shower and f^, the primary energy, la the 
typical energy emerging from the linear shower. The spatial 
dependence exhibits the large non-llnear growth of the pairs 
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possible at s/8q« 1+X"^. Dependins on the value of X, the srowth cay 
be ejcplootvc (X »1) or adiabatic (X«l). Thlo paraaater eatlcated 
for the galactic center would loply a alow buildup of pairs. An 
Interesting aside lo the posalble application of the explosive 
eolutlone to gaoaa ray bursts. In any case, the completed non-linear 
ehower has now Increased the particle flux frora ~1030 o±/a to ->10'*3 
e±/s at a distance of **10^5 cq. 


Electrons and positrons will now drift Into a region^ possibly 
an HII cloud'2^ and annihilate to fora the 0.511 HaV line. Palro at- 
now past the point whore the block hole could gravitationally 
redshlft the line out of the observed range and are also cool enough 
not to significantly contribute to the compact radio source past 
cn. It can now been seen that variability In the 0.511 line 
could be caused by variability In the dynano Itself. For Instance 
pinching of the beaa could cause the observed changes of line 
Intensity and would predict correlated hard X— ray and 0.511 KeV line 
variability. 

The density of pairs leaving the bean Is ~10'^ cra“*3. if nn HII 
cloud Is to act as a beam bag and subsequent annihilation region, 
the energy density of the beam should be comparable to the energy 
density of the cloud. If the pairs have energies of 1-10 MeV, a 
calculated cloud density Is "cloud cn“3. This Is 
consistent with H 2 C 0 { formaldehyde) and NH 3 Umnonla) observations. 

SUMMARY 


The model developed here supplies an annihilation region with 
the correct flux of relatively cool pairs at the proper distance 
frora the engine. Hard X-rays would be produced with comparable 
luminosity to the 0.511 MeV line and their variabilities would be 
directly correlated. An HII cloud appears to be consistent with 
stopping cooled pairs from the Jet. The proposed jet at the galactic 
center deposits more of Its energy In the 0.511 MeV line as opposed 
to the radio unlike Its more powerful counterpart In an active 
galaxy. This could be telling us the position of the acceleration 
point with respect to the central material. For an active galaxy 
acceleration occurs further out; the kinetic energy feeding a radio 
source. For an object like the galactic center, acceleration occurs 
closer In; the kinetic energy converting to rest mass which feeds an 
annihilation line. 

ACKNOWLEDGEMENTS 


The author would like to thank R. Lovelace, D. Lei ter, R. 
Ramaty, R. Llngenfelcer and the High Energy Astrophysics Theory 
Croup at Coddard for useful discussions. 


5 



ORIGINAL PAGE IS 
OF POOR QUALITY 


REFERENCES 

1. M. Levenchal and C, J. MacCallu.i>, The GalacClc Center, AIP 
Conf. Proc. 03,132(1982). 

2. K. Y. Lo, The Galactic Center, AIP Conf. Proc. 83,1(1982). 

3. J. H. Lacy, The Galactic Center, AIP Conf. Proc. 83,53(1982). 

A. J. L. Hatteaon, The Galactic Center, AIP Conf. Proc. 

83,109(1982). 

5. I. Gatley, The Galactic Center, AIP Conf. Proc. 83,25(1982). 

6. A. S. Jacobson, The Galactic Center, AIP Conf. Proc. 
83,123(1982). 

7. R. V. E. Lovelace, Nature, 262,649(1976). 

8. R. V. E. Lovelace, J. MacAuslan and M. Burns, Particle 
Acceleration Mechanisms in Astrophysics, AIP Conf. Proc. 
56,399(1979). 

9. Robert L. Brown, Ap. J. ,262,110(1982). 

10. R. D. Blandford, The Galactic Center, AIP Conf. Proc. 
83,177(1982). 

11. I. D. Novikov, these proceedings. 

12. R. E. Lingenfelter and R. Raiaaty, The Galactic Center, AIP 
Conf. Proc. 83,148(1982). 

13. M. L. Burns and R. V. E. Lovelace, Ap. J. ,262,87( 1982) . 

14. M. L. Burns In preparation. 

15. R. Gusten, The Galactic Center, AIP Conf. Proc. 83,9(1982). 


6 

JHJihTjOfjoai 


J 






ORIGINAL PAGE IS 5. N 8 3 27930 

OF POOR QUALITY 


COMPARISON OF PHOTON-PHOTON AND PHOTON-MAGNETIC FIELD 
PAIR PRODUCTION RATES 


M. L. Burns* and A. K. Harding 
NASA/ Goddard Space Plight Center, Grocnbelt MD 20771 


ABSTRACT 

Neutron stars have been proposed as the site of gaoma-ray burst 
activity and the copious supply of MeV photons admits the 
possibility of electron-positron pair production. If the neutron 
star magnetic field is sufficiently Intense (> lO^^ G), both photon- 
photon (2 y) and photon-magnetic field (ly) pair production should be 
important mechanisms. Rates for the two processes have been 
calculated using a Maxwellian distribution for the photons. The 
ratio of lY to 2y pair production rates has been obtained as a 
function of photon temperature and magnetic field strength. 

TNTnnnnnTiON 


Observations of the spectra of gamma -ray bursts indicate the 
presence of significant numbers of high energy photons the *iev 
range. Some spectra have features at energies between 350 and A50 
keV*, which have been Interpreted as red-shifted annihilation 
lines. The amount of ’ the shift (-10%) is the gravitational red- 
shift expected from the potential well of a neutron star. In 
addition, absorption features have been observed in many of the 
spectra in the region 20 - 60 keV, which if Interpreted as cyclotron 
absorption, indicate magnetic fields of order 10 G. The evidence 
seems to suggest that the emitting regions of these sources are nea 
the surfaces of strongly magnetized neutron stars. 

If this is the case, then one-photon as well as two-photon 
processes might be expected to contribute to the production of pairs 
in gamma-ray burst sources. One-photon pair production, a 
order process which is forbidden in free space, is allowed in the 
presence of a magnetic field. If the photons have a thermal 
distribution, then significant pair production rates will occur when 

(kT/mc2)(B/Bgp) > 0.1. 

ONE-PHOTON PAIR PRODUCTION RATE 


We first calculate the rate of lY pair production in a hot 
photon gas where the photons have a Maxwellian distribution. The 
rate for a single photon with energy E propagating at an angle 0 to 
a constant, homogeneous magnetic field of strength B Is J 


*Also University of Maryland, College Park, MD 207A2. 
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r 


YB 


_l CO BotnO 

2 a 


cr 


T (X). 


T(x) - 4.74 x“^^^ Al^ (X"^^^) 


(l> 

(la) 


~ 0.46 exp [■’4/3x1. (IB) 

- 0.60 X»1 (Ic) 

where x = (vez Kl — ) B ■ 4.414 x lO^^G, and Al lo Che Airy 

function. cr 

Theae are Che cayaptocic expreaalona in Che liaic where Che 
quanCua nuabera of Che cagnecic field pair otaCca are large. We 
will dioeuoa Che region of validity of Chia expreoaion below. The 
pair production rate for a diacribucion of photona. 0(E,T), will be 

Ryb(B.T) - 2» / ainOde / dE ♦ (E,T) r (E,B,9) a~^ (2) 

® \ln 

where ■ 2ac2/oln6 la the threohold photon energy. We take a 

Maxwellian dlotrlbution for the photons, nortaalized to constant 
photon nuaber density: 


«(E,T) - 1^3 exp (-E/kT) 


(3) 


If we isake the further approxicsatlons, x ^ ^ 2 qc^ 

then an analytic escpresslon can be obtained for the pair production 
rate. The E integration of Eqn (2) can be perforced by the laethod 
of steepest descents. The integrand has a saddle point at * mc^ 
[(8/3)(kT/cJC^)(B /B slnB)]^'^, about which the cajor contribution 
to the integral lo located. The 9 integral can then be perforoed by 
noting that the integrand peaks very sharply around slnO * 1. The 
result for the photon distribution of Eqn (3) is 

2 2 B 

R^b^B.T) a 7.92 X 10 ^ exp [-2 ^ ° W 


where the above approximations translate into the regions of 
validity for this expression: 



L 

B 


1 . 


cr 


kT . 3 B 

^ 2 B 

cr 


(5) 


It is evident that Eqn (4) is not valid for magnetic fields 
approaching the critical value or at low temperatures, where most of 
the photons have energies near threshold. 
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To obtain the behavior of the one-photon pair production rate 
in the full region of intcrcat, we have nuaerically integrated Eqn. 
(2), taking into account the photon energy threahold end the full 
aflycptotic rate [Eqno. (1), (la)J. Figure 1 ahowa thla rate per 
photon as a function of kX end B. Due to the csponentlal behavior 
of Eqn. (lb) at low photon energiea, the calculated rate la a very 
aenaltive function of T and B at low tenperaturca. At tcasperaturea 
below kT ■ nc^, alight varlaclona in B can change the rate by oaay 
ordera of eagnltude. 



Fig. 1. lY pair production rate per photon versus kT/mr.^ for 
different magnetic field strengths. 

RATIO OF ONE-PHOTON TO TWO-PHOTON PAIR PRODUCllOH 

The pair production rate via the two-photon process cr-n be 
calculated for the same photon gas with a Maxwellian distribution of 
energies In order to directly compare the one-photon to the two- 
photon rate. In the case of photon-photon pair production, we oist 
Integrate the cross-aectlon and the photon distribution over both 
photon energies, E and E* , and the angle, 0, between their 
propagation vectors (see eg. Ref. 3): 

n » « 

R (T) ■=■ 2nc / d0 sln0 (l-cos0) f dE f dE' 

. (B,T) . (E',T) a (E.E',0) a"^cn^ (6) 

2 4 

where E . (E>0) ■ 2 ra c /E(l-cos6) 

rain 
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The cross'-sectloa ie** 

o(E,E*,0) - I e(T) (7) 

vhere 

+ At - 8) ta - (t + 4) I 

/t - /t^ 

end t(E,E*, 6) - (Wooe). 

end the photon dlotrlbutlon le given by Eqn (3). IhtserlcaL 
Integration of Eqn (6) then enables us to evaluate the relative 
Itaportance of the one- and two-photon processes* Figure 2 ohosrs the 
ratio of the one-photon to the two-photon pair production rates an e 
function of kT. The vertical scale plots, on the left hand aide, 
the photon density (which Is kept constant with T) for vhlch the two 
rates are equal aud on the right hand side, the cctual ratio of the 
rates at a flsed density ol lO^S cn~3. Also plotted Is the 
blackbody photon density n^jCdashed line) which lo the tasicun 
photon density echievable at a given tcsperature. At tenperatures 
below ■ mc^. which are the tenperatures of Interest for gesna-ray 
burst sources, a change In the nagnetlc field of one order of 
mgnltude corresponds to taany orders of cagnltude In the ratio of 
the two rates* Frcs Figure 2. one can obtain, for a given photon 
density and nagnetlc field strength, the tcsperature at which either 
procees donlnates* For execple, for a photon density n - lO^S ca“3 
and B * 10^2 one-photon procees donlnates over the two-photon 

process at kT > tac^* For B > 4 x 10l2 one-photon pair production 
Is the dominant process at all temperatures* 

kT (keVl 


10* lO’ 10* 



Fig* 2* Ratio of ly to 2y pair production rates for Maxwellian 
photon distributions at teaperature kT* Curves are labeled with 
values of the magnetic field in units of 10^2 q . 


10 




Discossioa 




Pooh 




• vjr 


^^ALriy 


In a cagnctlc field, electron and positron eneroles 
perpendicular to the field are quantised, f^th the energy separation 
betveen the lerela increaalog^ vith field etreasth. Ixhea the photon 
energies are not large ccepared to the cpaclng betL^een these Landau 
states, the pair production rates are significantly Influenced by 
quantun effects* In the preceedlng calculations, ife h&ve neglected 
the effect of the discrete c4e- ototes on the pair production 
rates* In the case of the one-photon rate, quantua effects are 
Icportant uhen the quantity 2 (E/2nc^)V(B/B ) Is scsllS (le* for 
Iw temperatures cad high field strengths), uecrcaclng the rate 
belof the value given by the asymptotic expression [Eqa (lb)]* The 
tt#o-photon pair production rate has not been calculated In a 
cagnetlc field, so ue have need the free^pacc rate, vhere the 
electron and positron states are aoouaed to hz plane iraves* The 
effect of the cagnetlc field nl^t be eetlcated by exaolalcg the 
behavior of the Inverse process, tuo-photoa annihilation, which has 
been calculated In a strong cagnetlc field* ^ There is no 
significant deviation froa the free-nipace rate until the field 
approaches 10^ ^ C, and at this point, the one-photon annihilation 
rate begins to donlnate. Therefore, froa the proceeding arguesnt, 
our calculation of the ratio of the two rates is probably cost 
accurate below B 10^ ^ G and at higher temperatures* 

Photon distributions other than llaxwelllan cay be core 
realistic, since gacna-ray burst source csalctlng regions are not 
likely to be in equlllbrlua* The calculation presented here is ceant 
only to give an Idea of the relative behavior of the one- and two- 
photon processes* The actual rate of pair production In the source 
region would be an cqulllbriua solution Including other processes 
affecting the distrl^tlon of photons and pairs, such as 
annihilation, synchrotron radiation and absorption, and Coepton 
scattering* Cocparlsons such as this, however, cay aid In the 
construction of the rsorc self-consistent oodels by indicating which 
processes are important* 
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PAIR PRODUCTIG?! Af® ANNIHILATIOM IM STRONG HAGNETIC FIELDS 

3# K« Osughorty 

Unlverelty of North Corollna^Achovillo, AcheviXlo, NC 28804 

A* K« Harding 

NASA/Goddard Space Flight Center, Greenbolt, HO 20771 

ABSTRACT 

Electror*agnetic phenoisena occurring in the preacnco of otreng 
oagnotic fiolda cro currently of greot intereot in high^-enargy 
aotrophyoico* In particular, the process of pair ^^oduction by 
single photons in the presence of fields of order 10^^ Gcuos io of 
ieportsnee in cascade codeia of pulsar gsnea ray e:siooion, ond cay 
oloo becon^ significant in theories of other rediotion phenomena 
whose sources cay bo neutron stars (e.g*, gacca roy bursts)* In 
addition to pair production, the inverse process of poir 
ennihilotion is greatly affected by the presence of cuperstrong 
magnetic fields. The cost significant departures from annihilation 
processes in free spcco are a reduction in the total rote for 
annihilation into two photons, a broadening of the faniliar 511-keV 
line for ennihilotion at rest, and the possibility for annihilation 
into a si^^le photon (whi(^ docinatco the two-photon annihilQtion 
for 8 > 10 Gauss). The physics of these pair conversion processes, 
which is reviewed briefly, con become quite complex in the 
terogauss regime, end con involve colculations which ore technically 
difficult to incorporate into models of emission cschcnisms in 
neutron star magnetospheres. However, recent thcoreticuL work, 
especially in the case of pair annihilation, also ouggesta potential 
techniques for core direct ceasurecents of field strengths near the 
stellar surface. 

INTRODUCTION 

The observational discovery of pulsors in the late sixties was 
rapidly followed by their identification with rotating magnetic 
neutron stars. Early models of these objects, in %diich the magnetic 
fields were postulated as the means by which the otcllor rototional 
energy could be converted to electromagnetic rodiation, in feet led 
to enj^mous estimates of the surface field strengths, on the order 
of 10 Gauss or even more. Macroscopic fields of such intensity 
are roughly six orders of magnitude greater than the strongest 
attainable laboratory fields (which are currently the megagouss 
fields generated jy implosive flux-compression techniques^). In 
fact, the field strengths thought to exist in the interiors and 
magnetospheres of neutron stars nay well be the highest values 
occurring in nature. 

In light of the fact that teragauss magnetic fields ere so 
exotic by cerrestriol standards, it is not surprising that the 
theory of electromagnetic phenomena occurring in such on environment 
is for from complete. Although many of the fundamentals in this area 
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wore inveotigoted well before tha dlocovery of puloaro^^ oven these 
early Inveotigotiona hove only gradually bocooe well known to 
estrophyoiciote. Hence a considaroblo ccount of effort hao been 
apcjnt, oopeclolly in the last decode, on offorto both to Increooo 
^r undorstonding of quontum olectrodynaaico In strong fiolda, end 
to Incorporate this knowledge into opoclflc oodela of neutron otor 
nognetospherea and their cfniooicn rRochonians. 

In this diocuooion we will review coRie recent devoloptncnta in 
these orooa, especially thooe relating directly to eloctron-pooitron 
pair conversion processes. The treatinent is brief and incoriplete, 
but is intended to convey roughly the lovol of our current knowledge 
of these processes, and to indicate seise of the difficulties which 
this sort of physics presents to those interested in building nodela 
of pulsar cascade showers and other neutron star ccaission 
nechonisino. 


QUANTUM ELECTRODYNAMICS IN STRONG MAGNETIC FIELDS 


Although on a acole of the dinsensions of o neutron star the 
magnetic field is certainly nonuniform, the length scale of interest 
in quantum electrodynamics is so much emaller (on the order of the 
Compton wavelength) that in calculations of electron-photon 
interactions the field may bo considered perfectly uniform and 
infinite in extent. It may be noted that this is already a good 
approximation for typical eccolorator magneto, ond chould in fact be 
much better for neutron otor dimensions, whether theee fields might 
be simply dipolar or of some more complex multipolar form. Moreover, 
the field la believed to be so intense that its treatment os s 
classical or prcocrlbed "oxtornal'' field (which io not itoelf 
influenced by the particle interactions) is also en excellent 
approximation. Hence the fundamental tool needed for calculating 
electromagnetic processes in neutron otor mognetosphores ia the 
quantum mechanical solution for electron/pooitron motion in o 
constont, uniform magnetic field. Fortunately this la one of the 
coses for whic^ exact solutions of the reletivlstic Dirac equation 
are availeble . He will not discuss the details of these 
wovefunctions here, but will concentrate on the etiergetics and 
kinematics easocioted with them. 

The energy levels of a Dirac electron moving in a uniform 
magnetic field B may bo written in the form 


r 2 2 

s [c p 


m^c^d 


2n3/B^^)J 

CP * 


1/2 


( 1 ) 


where p denoted the component of, momentum parallel to the field 
® ^ 10^ Gauss, and n = 0, 1, 2, ... Tha 

form or this equation shows that the porollel momentum is not 
affected by the presence of the field, while the trsnsvorso motion 
is quantized. The critical field strength B is seen to be a 
combination of fundamental constants whoso vo%e is ouch that a 
tronsition between adjacent orbitals produces on energy change 
comperable to the rest maas of the electron* 

It is worth noting that the fully relotlviatlc Olroc solutions 
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to tho cqxioticns of eotion, which correspond to tho ensrgy level 
forculo given cbovot should be proforrod to tho corrooponding 
nonrolotiviotic SchrSdingor wcvofunctiono in soot cpplicationo to 
neutron ctor estrcphyoico. Thio io true in particular when either 
the psrollel c^ssntua p or tho product n3/Q^^ boco::^a oignificant 
coeporod to tho rest energy tsc*. For fiolSj for bolow 8 f tho 
trensverso energy bocon^s lorgs only for cnorcous values or*^n (in 
which caoe the trenoitiena to lower levelo produce what io usually 
tercied oynchrotron radiotion). Kowsvorf at field strengths 
conporcblo to B , oven nodcroto voluos of n (norr>ally ossocioted 
with nonrolotiviotic "cyclotron" radiation) can produce offoctivoly 
relativiotic behavior. 

Sons obvious chorectoriotica of the roletiviotic notion oay bo 
inferred froa Figure ly which shows tho firot few energy lovols in a 

field of 5 tcragsuso. (Horo 
tho porollol Gonantim p hoo 
boon cot to zero.) As the 
figure indic'atoo, eoch 
energy lovol above the 
ground state Eq octually 
corresponds to two distinct 
spin states, which cay be 
thought of GO "up" or 
"down" olong the field 
direction. It cay be found 
from tho relotivistic 
Dirac wavef unctions that 
transitions Involving spin 
flips ore leas proboblc 
then the corresponding 
"no-flip" tronaitions. 
However, the full set of 
spin states cust be 
considered in detorcining 
tho rclotive populations of 
excited states in neutron 
star magnetospheres. In 
addition, the decreasing 
spaclngs between levels for 
increasing quantum numbers 
n implies thot multi-level 
populations will produce 
not only multiple hormonlcs 
of emission lines, but also 
line-splitting effects at 
each harmonic. 


In terms of perturbation theory calculotions, pair production 
and annihilation effecta in strong magnetic fields are described not 
only by _the familiar second-order Feynman diagrams for these 
processes'^, but also by first-order diagrams involving only a single 


B = 5 * 10^^ G 


-^4 


-^3 


H2 


E-E« 

fXJH 


S=*l 


S=-l 


Fig. 1. Energy-level diogron for a 
Dirac electron in a uniform 
magnetic field of 5 teregsuos. Here 
n la the orbital quantum number as 
given in cq. (1), and the parallel 
momentum p s 0. Scale at right is 
the energy spacing in units of 
cyclotron energy hw,, = dhB/mc. 
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photon vertex. However* the Rsoning of "firot-order'* hero must be 
clorified. In contrcot to quontuni eloctrodyncmica in free spoce 
(where in fact ell firot-order procccoes oro kinenatically 
forbidden), tho Oiroc electron wavefunction itoelf fully deocribea 
the interaction with the conatent, uniform tnagnotlc field to ell 
ordera, while tranoitiona induced by redlation-field photona are 
deacribad by perturbation 
theory. The diatinction ia 

ahown in Figure 2, which 
comparea the second-order 
process of Bremaatrehlung in 
field-free apace with the 
"firot-order” process of 

magnetic Brcmastrahlung 

(aynchrotron radiation), in a 
uniform magnetic field. (Note 
that there io a second 

contribution to c^oe (a), in 
which the vertices are 

interchanged.) In (b) the 

electron interaction with the 
macroscopic B-flcld ia depicted Fig. 2. Comparison of oecond-order 
08 cn infinite number of diogrom (a) for Bremaatrehlung in 
zero-frequency photons, as free apace with "first— order 

opposed to the single "photon" diagram (b) for magnetic 
associated with the nicroocoplc Breaestrahlung (synchrotron 

Coulomb field of a nucleus in radiation) in external field, 
case (a). 

Of tho four pooolble ways to draw. Foynmon diagrono for 
first-order transitions, two represent Just the fcnilior proccaoes 
of synchrotron radiation snd ebaorption. Obviously these effects ore 
oboervoble for fiolda for below the critical value B--# Although it 
should be noted that the behavior of these efrccto in tho 
strong-field regina io quite different from that associated with 
interotellor or accelerator fiolda • However, the remaining two 
first-order procosoca, which are seen to be pair production and 
annihilation, ore easentiolly quantum-mochanical effects and as such 
become significant observable) only for fields approaching 

B . 

The kinematics of oil the first-order transitions ere 
determined by two equations, one for conservation of energy and one 
for conservation of porollel momentum: 

•heu s Ej + (2o) 

fitoco8<?/c s p + q (2b) 

where (E.,p) and (E. ,q) ore the total energy end parallel momenta of 
the poaltron and electron reapectively, and 0 la the angle between 
the photon wave vector ^ and the field direction. There la no 
equation for the conservation of tranaverae momsntum becauae the 
field Itoelf partlclpatea In the tranaverae momentum tranafer, and 
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olnco tho field llneo are eoouasd to ba infinitely "rigid** it io not 
poooiblo to dotorsalno tho trcnovorco Eoacntusa o«chcngo oufforcd by 
tho portlcloo* It lo porhopo worth noting, however, that tho 
rigidity of tho flold llneo lo oloo required to guarontoo tho 
conoorvotion of total energy ensong tho portloloo (oi^uotion 2 q) • In 
tho CQOO of noutron otor risgnotio fioldsi tho oocun^tion of 
Infinitely rigid flold llnoo io on extrcaoly good cpproxlnDtlont 
einco tho flold lo both cuparotrdng end end ”enchorod" (mochenlcelly 
cupported) by tho Rtooo of the otor Itoolf* 

PAIR PRODUCTION 

Tho procoDS of rtagnotlc pair production (tho convoralon of a olnglo 
photon Into on oloctron-pooitron polr In tho prooonco of end 
extern^ I nognotlc flold) hoo boon Invcotigatod by o nunibor of 
euthoro^^^, end the esocntiol rcsulto hove bocn known (olthougl-i in 
oona cooes svollcblo only co unpubliched doctoral thooes) olnco tho 
early fifties* Tho quontity of phyolcol Interest in dooling with 
this procces lo tho ottonuatlon coofficlcnt (invoroo of tho noon 
free poth) for o photon of epccified energy OJ end polorizetion 
vector £ » propogating at ooma angle & to tho uniform mognotlc 
flold (Honcoforth natural units, in which -h = c = 1, will be 

assumed*} The octuol 

colculotion of the ottonuatlon 
coefficient, booed on the 
first-order Feynman diegrom 
discussed above, may bo 
performed In tho Lorentz freme 
for which the photon motion io 
porpondlculcr to the field 
direction, with tho 

undorstanding that tho results 
(for unpolorized photons ot 
least) moy be gonerolizod to 
srbltrory directions of 
propagation by making Lorentz 
tranaformationa parellol to B. 
Wo will follow this approach, 
Fig. 3. Kinematically ollowod but we must remember that tho 
ototes for tho created eVe“ pair final step of transforming the 
with quantum numbers (J,k) . roculta must eventually bo 

performed in oppllcctlona whore 
k*B = 0* 

Before wo give tho results of Ihe dynotnicol calculstiono, It 
will bo unoful to consider the kinematic aspects of the pair 
production process. The final states of this process (ncmoly tho 
pair) may be labeled by two quantum numbers (J,k) describing the 
energy levels of tho electron end posllron respectively. As 
discussed in the preceding section, tho kinematicsl rcqulromants may 
be expressed as equations (2o,b). Now for a given photon energy end 
field strength B, these oquotions may be solved for the parallel 
momentum p! 



t 
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2 2 1 2 
P ■ p(j»)c) - * n [o' - 1 - (j + k) B' + (j-k) — ^] ^ ( 3 ) 

4u* 

where <u*s cu/2n and B* s 8 /B . The ooergy-mofflontum roetrlctlone ore 

equivalent to the rcquireiaciW that be a nonnogotive number } for 
given C 4 I and By any pair atotoa (J»k) which moot this requirement 
will be allowed. Ao ohown in Figuro 3, the region of allowed atatoa 
in the (j,k) plane la thus encloaod within the line corresponding to 
the equation p(J,k|W,B) = 0. The aum of the individual tronaition 
rates for these final states will yield the attenuation coefficient 
for the pair converoion proceoa. In addition, the relative 
probabilities for the various transitions over the (J,k) plane 
determine the energy distribution of the emergent pair, which is 
also of importance in applications (e.g«, pulaor cascade modela). 

The explicit calculation of the attenuation coefficients in 
terms of the Dirac wavefunctions has been fully dlecusaed in the 
references mentioned above and will not be repeated here. The 
results for the cases of photons whose electric vectors ore 
polarized parallel and perpendicular to the field direction 
respectively are given by the expressions 

R — I Z -i- |(E E 4o*-p^H|H(j,k)| + |M (j-l,k-l)l ) 

• 25 j klp^j^l 5 k 


2 . 

+ 2/jk" B'n Ill*‘(j,k) M(j-1,k-l) + M(j,k)]} (4a) 


R m^l T ^{(ee -ha^^XlMCj-lrk) |*+|H(j,k-l)|^) 

^ j k'**jk‘ j k 


2 « 

-2/jk B'n tM^{j-l,k)H(p,k-l) + M^(J,k-l)H(j-l,k)]} (4b) 


where K “ • M(J,k) - (-i)' 

B* 


G-S fsi 


G-S 


G-S 


Cl 




and G - laax (j,k), S - nln (J,k) 
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At ripot glonce those oquetiona are not vary illtfitinating, but fpoa 
then we can irssodiotoly Infer ona irportont eopoot of the behavior 
of the ottcnuotion cosfflclento sa functlono of photon energy ond 
field strength. The 
eppeorence of the 
noraontun term p In 
the donoainator of 
each oui»sand irspllee 
that If ony p 
vanieheo, the entire 
exproaoion becon»a 
singular. Now from 
equotion (3) it may 
be seen that fof 
each integer pair 
(J>k)y only certain 
combinationa of 
and 0 can make p 
vanioh. This in turn 
iropliea that for 
fixed B the 
attenuation 
coefficients R(<s» ,B) 
are singular at a 
discrete sequence of 
energies , with 
each singularity 
resulting from 

particular (j,k) volues for which p = 0. In terms of Figure 3, the 
attenuation coefficients ere singular whenever co and B ere such 
that the lino p(J,k; co ,B) = 0 Intersects an integer pair (Jtk). 
Hence a plot of R{ co ,0) vs. co shows the sort of sawtooth behavior 
which is depicted in Figure 4. The energies at which the 
eingularities occur ore readily found from equation (3), and by 
plotting these energies it is found that the average spacing between 
peaks rapidly becomes smaller as the field strength is decreased, 
ond that for fixed B the peaks in successive fixed-length energy 
intervals dco become more numerous. 

This complex behavior was gotcd even by the earliest authors 
who investigated this process but for the maximum field 

strengths then considered to be of conceivable practical interest, 
the density of singularities in measurable energy intervals would be 
so great that only smoothed-out, average values token over each 
interval were considered to be of physicol interest. Hence 
asymptotic expressions for the attenuation coefficients, volid in 
the limiting regimes U) >> 2m and B « B , were derived from 
expressions (4). The cruciol steps involved *in this derivation are 
the replacement of the (j,k) summations by integrations over 
suitably chosen continuous variables, and the ^eterminction of 
appropriate^ j0symptotic forms for the generalized Laguerre 
polynomials ' . (For a recent discussion of this derivation, see 

reference'*' •) The final results are expressible in the relatively 



w (UcV) 

Fig. 4. Exact attenuation coefficients for 

photons propagating ot right angles to the 

field direction ond with polcrizotions 

parallel ond perpendicular to the field, 

plotted vs. energy for fixed B (here chosen 

to be B ). 
cr 
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, X = « 1 


( 5 ) 


A comparison of the asymptotic results with tha exact forms (4) for 
two sample field strengths in shown in Figure 5. 

In most aatrophysical applications to date^^c£gecially in 
models of cascade showers in pulsar magnetospheres * , the lost 

form of the attenuation 
coefficient has been used. 
Since the typical photon 
energies produced by 
curvature or synchrotron 
radiation in a cascade 
sequence are well above the 
MeV range t and the ambient 
field strengths ore usually 
assumed to be no more than a 
few teragauss, the complex 
behavior for near-threshold 
energies and fields B < B^^ 

has essentially been 
ignored. Moreover, the 
energy distribution of the 
pair has ust >lly been 
assumed to be gi/en simply 
by Ej. = CO/2. (It should 
also be noted that only 
unpolarlzcd attenuation 
coefficients are used in 
these models.) 

However, it turns out that this asymptotic form is often not a 
good approximation at all, since (as mentioned above) it is really 
necessary to Lorentz-transform the attenuation coefficients before 
they may be applied in frames for which ]<*B, does not vanish. Now it 
turns out that the transformetion low of the attenuation 
coefficients may be expressed as 

R (U/ B) > sin 6 R (u sin 6 . B) (6) 
— o “ 

This law expreaaes the fact that energetic photons, propagating at 
an angle & ^ 1/Ct) to the field direction, may be 
Lorentz-tronaformed down to energies near (or even below) threshold 
in the ” transverse” frame for which = 0* But it is Just such 
angles of propagation which are typical for the curvature-radiation 
photons which should inititote the pulsar cascade process. Hence a 
proper treatment of pair production by such photons must take into 
account the neor-threshold effects, and the resulting modifications 
to the estimated multiplicities (ratios of secondary pairs to the 
number of radiating primary electrons drown from the stellar 



Fig. 5. Comparisons of exact 
(unpolarlzcd) photon attenuation 
coefficients vs. asymptotic forms 
for 1- and S-teragaucs fields. 
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surfaco) nay have o olgnlflcant effect on the entire oaecedo 
development* Moreover, it turna out that tho energy diatributiona of 
the prtira in tho noar-thrcahold or high-fiold rcginsa, oa doteminad 
by the individual contributione to eKpresaions (<%), can broaden 
conoiderobly end even chow double-pasked behavior oa ehotfn in Figure 
6. (The double peaka imply that one meebor of the pair tends to got 
most of tho photon energy, while the other emerges much clocor to 
ito ground ctote.) Konce tho usual escur^itlon that E* a CO/2 io oloo 
likely to bo o poor epproxlnatlon. This io unfortunately only one 
example whore noro careful treatmanto of the underlying proceasea in 
neutron otor magnetosphoroa may force oignif leant revisions to 
current modale of their emission mechaniems. 


OWGINAL 
OF POOR 


Fig. 6. Energy distributions for 
one member of the crested pair, 
computed from an integration over 
tho (J>k) probability 
distributions. The diotributiona 
shown are oil normalized to unity, 
and ore plotted vs. pair energy 
divided Jjy photon energy. 


'0 0^ I 
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PAIR ANNIHILATION 

The inverse process of first-order, pMr production, namely poir 
annihilation into a alngle photon ’ , ia also of interest in 
neutron star astrophysics, but in this case the process becomes 
significant only for fields very close to B .In weaker fields the 
first-order process is dominated by the an^’log of annihilation in 
free space, in which two or more photons make up the final state. 
The two-photon annihilation process, which will be discussed later 
below, is itoelf strongly modified by ambient fielda of a few 
teragauas or more. 

The kinematics for one-photon annihilation are similar to 
equations (2) for poir production, where the roles of initial and 
final states are now rev»*rsed. These equations imply that for 
annihilation by pairs at rest, the emergent radiation must form a 
flat fan beam at right angles to the field direction. For thermally 



21 






ORIGINAL PAGE fS 
OF POOR QUALITY 


or otherwise broadened eVe" parallel mosnentuRi distributions, the 
fan besn would fill out to a form such os that shown in Figure 7. 
Corresponding to the angular distribution, the opoctral distribution 
of the photon, assuming emission at specific angles, have 
characteristic asymmetrical shapes as shown in Figure 8. Both the 
overall shift toward higher energies and the increased broadening 
are due essentially to the Doppler effect caused by annihilation of 
pairs with nonzero net parallel momentum. 



ti«(UeVI 


Fig. 7. Angular distribution of Fig. 8. Differential one-photon 
the one-photon annihilation annihilation spectrum emitted 
radiation resulting from from a ground-state gas of 
Gaussian electron-positron electrons and positrons with 
distributions with momentum Gaussian parallel momentum 
widths of 50 keV, in a field of distributions, here with widths 
10 terogauss (see reference 17). 'Jf 1 (reference 17). 

Separate spectra are shown for 
three angles of e/aission. 

The kinematics for two-photon annihilation ore laore complex but 
are still characterized by a loos of transverse momentum 
conservation. As in the cose of the first-order processes, the 
electron wavefunctions (and here the propagator as well) correspond 
to the Dirac constant-field wavefunctions. Several dynamical and 
kinematical aspects of two-photon annihilation turn out to be of 
special interest for neutron star astrophysics. The first is a 
reduction in the total annihilation rote os conpared with the 
free-spoce value, accompanied by a flattening of the isotropic 
angular distribution of the emitted photop^ (tending toward a fan 
beam perpendicular to the field direction) • In addition, the role 
of the field as a transverse momentum absorber leads to a p^arply 
field-dependent broadening of the annihilation spectrum . In 
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h at reot, the cpectrua io no longor 

conotrainod by wmontum conaorvation to the foralHor 511-keV liL, 

aonaitlve both to field 
IjJ^hrn T**!**"^^!® omission (see for example Figuroa 9,10). The 

obaLvAhl°e®ri"rf ®ep®®i®lly intcroating in that it might bocome 
wervable to detectors with high energy resolution for fields of 

^”®9®®®®‘ Hence as a potential means for direct 
^n!hn::??nn ®trengths, the line-broadened ?i4hot®n 

• • h have wider usefulness. In this context 

.. should be also be noted that the angular dependence of the 

should, in beemc emerging from local "hot spots" 



Fig. 9. Two-photon annihilation 
spectrum for pairs at rest, 
observed perpendicular to the 
field direction, for several 
field strengths above one 
teragauss. The sharply 
field-dependent line broadening 
is evident and is here due 
entirely to the loss of 
transverse momentum conservation 
In the intense field (i.e., no 
thermal broadening is included) . 



Fig. 10. Two-photon spectrum for 
annihilation at rest in a field 
of 10 teragauss, as viewed from 
several polar angles. 


F^inally, it should be noted that the range of field strenoths 
o..r wh.ch the line broadening end overall ?.te deirltf 0 ^!;. 
«m PHoton annihilation process becomes observable is really quite 
small, on the order of only one teragauss. The same statement holds 
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for the **croosover* field strength (cbcsut 10 teragauss) at »<hich the 
total one-photon cnnihlleticn rate begins to dominate Uw two-photon 
rate« as shown in Figure 11« Given this eKtreme sensitivity to the 
ofiibient field strength, it may eventually becone possible for 
detectors with both hi^ spectral end temporal resolution to provide 
accurate locol values or even siaps of the surface fields through the 
measurement of cnnihilotion radiation* 



fig. 11- Comporison of the total 
rates for annihi lotion into one 
end two photons at rest, as a 
function of field strength- The 
dashed line corresponds to the 
limit for two-photon 
annihilation at rest in free 
space, taking into account the 
fact that the electrons ore in 
this case polarized as they 
would be in the ground state if 
the field is present- 


CONaUSION 

The detailed physics of processes in tcragauss magnetic fields, 
including the sicsple first- and second-order processes of the type 
discussed above, is still far from complete. The calculations of 
even the first-order processes involve a great deal of technical 
difficulties which ore especially troublesome in the conteKt of 
neutron star astrophysics, where often time-consuming calculations 
are already required to estimate emission spectra from the models- 
Indeed, it appears in many cases impractical to incorporate "exact" 
expressions such as equation (A) into models ot pulsar cascades or 
garnma-ray burst eources, ot least with computing resources that are 
less than what curroitly fall into the cate^ry of supercomputers- 
On the other hand, the analytical difficulties ond pitfalls 
associated with attempts to improve or extend the asymptotic results 
for these processes are equally formidable, and the reliability of 

the results is open to question* 

These seme problems ere only aggravated if we begin to consider 
various other fundamental processes, most of them intrinsically more 
complex than the sim4)le pair conversion effects discussed here, 
which may also play significant roles in neutron star astrophysics. 
Among these may be, pmentioned such eff^-cLg as the field-induced 
index of refraction^®*-^^, photon splitting , trident production , 
and such complications to the pair annihilation process as 
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••nagnetlc* positrcnium forcatim end ennihilotion of posltrcno ttLth 
the otafilc elcctrone in the eurfcce loycre of the neutron oter* At 
preo^^tf it io safe to oey that the theoretical probl^us in this 
orea Mill only olowly be put to rest. 

At the GGce tice, it io encoureging to note that cven^ fairly 
clesjsntory taodelo of such proceeees oo pulsar cascade choMoro^*^ have 
already yielded reasonable agreeaent with the currently oveilcble 
observational data, in spite of the fact that the fundancntol 
cenvereion processes have been treated in rether crude wayo to date. 
Hence it Cisy be that refine^icnto to the theoretical treatcento of 
these processes Mill only intprove this ogreenjent end not produce 
drastic or qualitative changes in the existing c^dsls. At present, 
hoMever, it appears necessary not only to ctakc the usual pica for 
core observational data, but also to csphacize the need for further 
theoretical investigations of the elcssentary processes of high-field 
physics end their careful trcotcient in models of pulsars, pulsating 
X-ray sources, end gsnnia-roy bursters. 
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ABSTRACT 

For a satisfactory understanding of astrophysical annihilation 
radiation, cxpeclally that observed froa the galactic center 
direction, the Interaction of positrons with the ambient medlua 
roust be carefully Investigated. Although hot, Ionized regions nay 
be Important sources of annihilation radiation. In this report I 
will exanine nalnly the simpler processes occurring In 
low- temperature neutral hydrogen gas. The goal Is to set Holts on 
conditions In the annihilation region by using the predictions of 
atonic theory compared with the observed y-ray line width, 
continuum strength and time dependence. 

IHTRODUCTIOM 

Suppose that positrons of energy high compared to the 
Ionization energy of atomic hydrogen (13.6 eV) enter a region 
containing cold, non-lonized atomic hydrogen and nothing else. 
Assume further that no magnetic field or radiation Is present In 
the region and that the hydrogen density Is much less than 10^2 
cm”3. He ar« thus considering a very simple. Idealized model for 
the galactic center annihilation source and will try to predict the 
properties of the annihilation radiation emerging from this 
source. In spite of the simplicity of the model, it will be seen 
that a good deal of interesting physics is Involved in its 
analysis. At the present time, there are no experimental data on 
the collision of positrons with atomic hydrogen, as there are for 
molecular hydrogen* and atomic helium. 2 On the other hand, the 
theoretical situation, although still Imperfect, Is relatively 
good, as I will show. 

The observations2-5 of the galactic center annihilation 
radiation are easily summarized: 

1. A narrow line is seen; the best measurements of its width 
(FHHM) gives 3.13 * 0.57 keV. With a detector resolution of 2.72 
keV, the source line width is 1.6(+0.9,-1.6) keV, consistent with 
zero. 

2. The energy of the line is very accurately that expected for 
e* - e“ annihilation at rest, 510.90 * 0.25 keV compared with 
511.00 keV in the laboratory. 

3. There is some evidence^ for the existence of a continuum 
component on the low-energy side of the line in addition to a 
power-law background. 

4. The intensity of the line radiation is definitely 
time-variable,S with a time constant less than about 1/2 year. 

In the next section the physical processes leading to 
annihilation of positrons will be described, and the resulting line 
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vddth wfll be shown to be too large to agree with observations. 

The standard codification, allowing the annihilation region to be 
partly Ionized will also be discussed. In the final section a 
possible loophole will be described, perhaps allowing retention of 
the original picture of a cold, non-lonlzed region and predicting 
sense observational consequences. Kuch of the discussion below Is 
based on the work of Bussard et al .^ and Crannell et al. ^ 

POSITROM ANNIHILATION 

Five different reactions can occur In positron-hydrogen 
collisions at codcrate energies (13.6 eV < E < 1 keV): 


♦ IJ 

(Elastic Scattering) 

(la) 

♦ H* 

(Direct Annihilation) 

(lb) 

♦ H* 

(Excitation) 

(ic) 

♦ e- ♦ H* 

(Ionization) 

(Id) 

♦ H* 

(Positronlum Formation) 

(le) 


SteckerS pointed out very early that In this energy range only 
the last process Is effective In annihilating positrons, since Ps 
forcation has an atosslc-slzed cross-section of order agE (-I0“16cc2) 
while process lb Is radiative, of order rg2 (-io-26cq2). 

Essential to this argucent Is the low density of scatterers; once a 
Ps atom Is foraed It has no further collisions In the short time 
before It annihilates into 2 v-rays (1.25 x 10~10 sec) or 3 
Y-rays (1.4 x lO”^ sec). (The opposite occurs in most laboratory 
experiments where high densities are maintained; Ps formation In 
this energy range Is merely an Inelastic process, since the Ps atom 
Is almost always colllslonally ionized before It annihilates.) 

With these considerations In mind, I can describe the life 
story of a positron very simply. It enters the cold atomic 
hydrogen region at son» high energy and loses energy by processes 
(Ic) and (Id). At about 100 eV process (le) begins to become 
non-negllglble, and for the rest of Its lifetime every collision 
carries a certain Increasing risk of Ps formation and Immediate 
annihilation. It is the velocity distribution of the Ps atoms at 
the Instant of their annihilation into 2 t-rays that determines the 
width of the observed annihilation line. Before examining this 
question more quantitatively let us first note that positrons of 
energy E produce posltronlum atoms of energy E - Eg (where 
Eq = 6.8 eV, the threshold for process (le) giving a rectangular 
distribution in the line-of-sioht component of veloctty v,. 


u t > wi I III uiitc II I i un u v. u u i vc i ul. i cjr v ^ , 

ranging between v- « * ^(E-Egl/rngl Since the Doppler shift 
of one of the annihilation y-rays is a = (v2/c)mgc2, one gets 
a rectangular line profile of width 


r X 2a = 1.430^E(ev)-6.80 keV. 

For this monoenergetic positron distribution to satisfy the 
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experimental constraint r<. 2.5 keV Eq. (2) requires E < 9.86 eV; we 
will see shortly that thls^ls an unreasonably low energy. 

In Ref. 6, a' detailed Monte Carlo simulation Is carried out, 

In which an ensemble of positrons Is follcv;ed downward In energy 
until annihilation takes place. The resulting T-ray line histogram 
corresponds to a width of 6.5 keV, In clear disagreement with 
observation. Rather than repeating the details here, I will use a 
simple continuous slowlng-down model to describe the thennallzatlon 
of the positrons. First, however. It is necessary to review tlie 
status of positron-hydrogen scattering theory, beginning with the 
Ps-formatlon cross-section. 

This cross-section Is difficult to compute accurately. In part 
because of the unsymmetrical relationship between Initial and final 
states In Eq. (le); the center of mass coordinate of the Ps atom Is 
not a natural one for the Initial e^ - H state, nevertheless, a 
number of two-channel calculations have been carried out both for 
s-waves^ and p-waves,*^ and they agree in predicting a great 
reduction below the Born cross-section** for Ps formation, but 
they do not extend up to energies of Interest here. For 'that 
reason, a phenomenological extension to higher energies was carried 
out,^*^ which used a modified Born approximation whose LaO and Lsl 
partial cross-sections were reduced to agree with Refs. 9 and 10, 
ano whose L > I terms were unchanged. The only real test of this 
approximation comas from a comparison with the total Inelastic 
cross-section of Hinick and Reinhardt, obtained by a 
sophisticated analytic technique. Although no distinction is made 
between Ps formation, excitation and Ionization, only the first of 
these Is energetically allowed for energies betweeii 6.8 eV and 
10.2 eV. There is reasonably good agreement In this energy range 
between the results of Refs. 12 and 13, encouraging me to use the 
results of Ref. 12 in the rest of the analysis. (Note that Ps 
formation in excited states has been neglected here; an increase of 
less than 20 percent might be expected in the cross-section at 
higher energy. ) 

I do not know of any positron-hydrogen ionization calculations 
(except for the Born approximation which does not distinguish e^ 
from e-.) A very simple analytic form has been devised by 
Lotzl** for the e“ - H ionization cross-section: 


01 _ 2.47 InE [1 - O.C e-0-56(E-D] 

E 

where E is the energy in Rydbergs. I will use Eq. (3) in the 
analysis, although there is no estimate of error in the e* case. 

On the other hand there is a recent close-coupling calculation 
of the positron impact cross-section for excitation of the n=2 
levels of hydrogen. up to at least E = 7 Ry («100 eV) these 
are well fitted by a formula like that of Eq. 3: 


®exc 

«ag 


2.73 (In (4/3 E) (1-0.63 e-0-531E) 
E 


( 4 ) 
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The n»2 excitation should dominate the total excitation process, so 
I will use Eq. (4) in the further analysis. 

In Fig. 1 the Ps-fomiation cross-section and the inelastic 
scattering cross-section (oj + opxc) plotted as functions 
of energy. To give an idea of the uncertainty in these quantities. 



Fig. 1. Positron-hydrogen cross-sections as functions of energy. 
The solid line is for Ps formation (Ref. 12) and the dashed line is 
oT ♦ oexc 3 and 4). The dots and the cross are estimates 

Of oj ♦ ogxQ obtained by cwibining the results of Ref. 12 with 
those of Ref. 13 and Ref. 16, respectively. 


I have added a few additional points: the six low-energy points 

represent the difference bewteen total inelastic from Ref. 13 and 
Ps-formation from Ref. 12; the point at 100 eV is a similar result 
from an eikonal approximation.^® These are both nigher than the 
cross-section used here, but not, I think, in too serious 
disagreement; for atomic research more accuracy is desirable, but 
for astrophysics the present results are quite satisfactory. 

Before carrying out a more detailed calculation of positron 
slowing down and Y-ray line shape I can draw some semi -quanti tati ve 
conclusions from Figure 1 itself. Notice that the two competing 
processes, inelastic scattering and Ps formation, are equally 
probable at a positron energy of 28 eV. Of the positrons surviving 
to reach this energy one-half will form Ps at their next collision, 
and from Eq. (2) they will give a line width of 6.58 keV. It is 
thus unlikely that the observed width of p £ 2.5 keV can be 
achieved. 

To account correctly for the large energy losses occurring at 
each inelastic collision a Honte Carlo simulation is needed.® 

But an approximate, qualitatively correct treatment involving a 
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continuous s1ow1n9*dovm opproxlaatlon is easy to fonnulate and 
yields results much like the more exact ones. Assume that all 
positrons enter with the sar.e high Initial energy Ei and lose 
energy according to the equation 


“ -Coexc (E) AE ♦ oi(E) All Nv(E) (5) 

dt 

where the cross-sections are from Eqs. 3 and 4, N Is the nu mber 
density of hydrogen and v Is the positron velocity, /ZE/nJa . The 
energy losses will be taken» as ae • 10.2 eV and at « 17 eV. 

At the same time, the swarm of positrons Is being depleted by 
Ps-forroatlon following the equation 



-ops(E) n(E) v(E) N 


( 6 ) 


These coupled equations can be Integrated to give the time history 
of a positron swarm or the time dependence can be eliminated: 


dn(E) ^ dn(E) dE 

* “3E 3T ’ 


(7a) 


« g(E) n(E), (7b) 

where g(E) = ops(E)/[oexc(E) ae ♦ oi(E) ai], Eq. (7b) 

|ias the trivial solution 

n(£) = n(E() exp f dE' g(E') , (8) 

J El 

and I have plotted n(E)/n(Ei) in Fig. 2, where Ei » 20Ry « 272eV 
is a high enough energy to be considered asymptotic, since g(Ei) 

Is very small. It is, of course, unrealistic to carry the solution 
below E« Aj, since there the continuous slowing-down 
approximation is grossly incorrect. 

The principal conclusion to be drawn from Fig. 2 is that the 
half-value energy, where only 1/2 of the original positrons still 
survive is at a high energy, E 1/2 = 39 eV. This Is a further 
indication that most of the Ps atoms formed are moving too fast to 
give the required narrow annihilation line. Furthermore, the exact 
line profile corresponding to this form of n{E) can be derived by 
integrating the rectangular line shapes of Eq. (2) normalized to an 
area proportional to dn(E). That is, 

P(A)oC dE' g(E')n(E') /Jr^E^ . (9) 

oEo^A^/mgC^ ^ 
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EMERGY (eV) 

Fig. 2. Percent of positrons originally at 272 eV surviving to 
reach a given lower energy while slowing down in neutral atomic 
hydrogen. — — 


^ as 


0 1 1 3 'V 5 ^ 

A (keV) 

Fig. 3. Gamna-ray line intensity as a function of distance from 
the line center a. The line is symmetric about a = 0. 
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In Fig. 3 I have plotted the annihilation line profile against a In 
keV. (The central parts of the line Involving very low values of 
positron energy, where the present approximation Is poor, have been 
extrapolated quadratically.) The width of the line calculated this 
way Is 6.75 keV, in good agreement with the Konte Carlo results of 
Ref. 6, and once and for all Inconsistent with the observations. 

The solution to the d11en?na Is usually presented as 
follows”: If one allows some 5-10 percent Ionization In the gas 

the slowing down of the positrons is so efficient that they do not 
form an appreciable amount of positronlum before they take up a 
thermal velocity distribution. In Fig. 4 a diagram from Ref. 6 is 
reproduced. It shows the thermal -average rates for the four 



Fig. 4. Rates {per unit target density) at which thermal positrons 
form Ps by charge exchange with H (Rce^^H^ O'” radiative 
recombination with free electrons (R^p/Og), and annihilate 
directly with free electrons (Rja/ne) or with bound electrons 
{Rda/nn), as functions of temperature. (From Ref. 6) 
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processes occurring in partially ionized hydrogen. At temperatures 
belov/ about 5 x 10^ K the dominant process is radiative 
recombination with an electron to form Ps followed by 
annihilation. The shape of the resulting line is 

P (a) = e-A^ATniec2 ^ dO) 

its width is^ 1.105 (T4)l/2 i^gv where T4 is the temperature 
in lO^K, and a satisfactory line width can be obtained. At 
higher temperatures charge exchange (Ps formation) with neutral 
hydrogen dominates; it rises very rapidly since only the tail of 
the thermal distribution above Eg = 6.8 eV can contribute. 

Once convinced that annihilation from a thermal sv/arm is 
occurring, we can, in principle determine the ionization fraction 
and the temperature if we know the line width and the Ps fraction. 

At present, of course, we have only an upper limit on the line 
width, while there is only rough and discordant data on the Ps 
fraction, which is obtained by a delicate process of curve fitting 
in the presence of a large y-ray background. In Ref. 4 the most 
probable Ps fraction was given as 92 percent (but consistent with 
zero). At this Workshop, however, Rieglerl^ suggested a value of 
about 20 * 20 percent. Clearly no conclusion is yet possible. If 
the lower value proves to be, correct, hov;ever, severe constiaints 
will be placed on the scenas^tST described above. In particular, it 
is not possible to keep the Ps fraction below about 60 percent 
while keeping the ionization fraction >5 percent. Several possible 
ways out of this dilemma would be to allow the annihilation to 
occur at high density^^ (to pick off the triplet Ps) or to flood 
the annihilation region with ultraviolet radiation (to ionize the 
triplet Ps before decay.) I will propose another possibility in 
the next section, by going back to the completely neutral case. 

A TIME-DEPENDENT SOLUTION 

In discussing Fig. 2, I did not mention the fact that about 5 
percent of the original high-energy positrons avoid being 
annihilated as Ps during the slowing down process. This conclusion 
was originally drawn from the Monte Carlo calculation of Ref. 6, 
but a similar estimate can be obtained from Fig. 2. Note that the 
continuous si owing-down model should not be carried to an energy 
below 23.8 eV, since below that point a single ionization event 
will bring the positron below the Ps-formation threshold. But 19.0 
percent of all the positrons reach 23.8 eV, and the probability 
that the next collision is an ionization event is Pi(23.8) = 
oi(23.8) / [oi * agxc * ®Ps^ = 0.1335. Similarly, the 
probability that the next collision is an excitation event is 
0.2776, lowering the positron energy to 13.6 eV. At 13.6 eV, the 
excitation probability is 0.1303; in effect we are doing a very 
simple Monte Carlo calculation for the last few collisions before 
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the Ps threshold is reached. The result Is 

n(E<5.8) - n{23.8) [Pi(23.8) ♦ Pexc^23.8) x Pexc<13.6)] 
s 3.2 percent. 


in qualitative agreement with Ref. 6. Three questions may be asked 
about this positron residue: how wide is the yny line it 
produces, what is the resulting Ps fraction, and what is its time 
dependence 

In this region {E < 6.8 eV) where energy loss is by elastic 
scattering only (&E/E » 2mg/mH per collision) and the 
annihilation rate is very Tow (due to collisions with H atoms), the 
continuous approximation is reliable. Using the almost exact 
e*-H annihilation and momentum- transfer cross-sectionsl« in the 
following equations 

dE _ -2mg CI^^T(E) Ev{E)N 
dt rnii 


= -aA(E) n(E) v(E)N (12) 

“3t” 

it is very easy to show that almost all the positrons cool rapidly 
and annihilate at very low evergy. (For example, 94 percent of all 
the positrons initially at 6.8 eV reach an energy of 0.05 eV, 
corresponding to about T = 600 K.) They will then annihilate with 
the bound electrons in the atomic hydrogen ground state; the 
quantum-nrachanical momentum distribution is the cause of the 
annihilation line width, since the energy of the positrons is 
nearly zero. 

If one has an accurate e -H zero-energy scattering wave 
functionTj"(^, r), where "Sf, r* are the co-ordinates of incident e 
and atomic e“ (in atomic units) respectively, the annihilation 
line profile is given by’^ 



r*” 



P(a)c< 

\ dq q 

I \d? jo(qx) (^> ^)| 

(13) 


J2A 



duriQC^ 

Humberston20 has carried out such a calculation, although with a 
different situation in mind, and from it an annihilation width of 
1.3 keV can be derived. (Several references6,21,22 have 
misquoted this as 2.6 keV.) This positron component gives a narrow 
annihilation line with a vanishingly small Ps component. (If some 
Ps is observed, a very small ionized fraction will account for it 
easily; at T = lO^K, ng/nH = 0.0025 gives 20 percent Ps 
fraction.) 
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In Ref. 6 we assumed that the broad line due to Ps formation 
and this narrow component were properly weighted and superposed, 
since we assumed a time-independent source. This coinbined 
annihilation line is not narrow enough to agree with the 
observations. The new suggestion I would like to make is based on 
the observed time dependence of the radiation. 

Suppose a short burst of positrons is injected somehow into 
our cold atomic hydrogen region. He would observe two successive 
annihilation phases: The positrons would slow down in a relatively 

short time and, as described, would form Ps in flight giving a wide 
line at first. Soon, however, the only remaining positrons would 
be in the component below 6.8 eV, annihilating slowly with a narrow 
line width. From Eq. (12) it follows that if this second component 
is being observed N > 10® cm"3 for a mean life of 1/2 year. At 
this density the first phase should take less than 10 percent as 
long, 23 about one or two weeks. The conclusion is that randomly 
timed observations of the galactic center are much more likely to 
see narrow lines than wide ones, and this may be the explanation of 
the present observational situation. 

It is important, then, to try to observe the galactic center 
source on a continuous basis, for the general purpose of charting 
its time dependence and specifically to look for the unique time 
dependence of the line width discussed above. 
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ABSTRACT 

A high^csolution observation of the active nucleus galaxy 
Ccntaurus A (NGC 5128) was made by the CSFC Low Energy Gasna-ray 
Spectrometer (LEGS) during a balloon flight on 1981 Hovcabcr 19* 

The taeasured cpectrua between 70 and 500 keV ic wll represented by 
a power law of the fora 1.05 x 10*^ (E/100 ph co*^ 

with no breaks or line features observed. The 98Z confidence (2o) 
flux upper limit for a narrow f< 3 keV} 511-kcV positron 
annihilation line is 9.9 x 10""^ ph s'" . Using this upper 

limit, the ratio of the narrow-line annihilation radiation 
luminosity to the integral > 511 keV luminosity is estimated to 
be < 0.09 (2o upper limit). This is compared with the measured value 
for our galactic center in the Fall of 1979 of 0.10-0.13, indicating 
a difference in the emission regions in the nuclei of the two 
galaxies. 


INTRODUCTION 

Centaurus A (KCC 5128) is a nearby galaxy 5 Mpe) with an 
active nucleus that is an Intense source of X-rays and gaenna-rays. 
The nuclear source is variable on tltae scales from days to years and 
is spatially unresolved at X-ray energies; the upper limit to the 
sire of the nuclear component of the X-ray emission is 0.3 arcsec^ 

(8 pc). Einstein observations indicate that above — 2 keV the 
nuclear component dominates the X-ray emission from the galaxy^. 

Th^re are three previously reported observations of Cen A by 
instruments capable of measuring a positron annihilation feature in 
the spectrum: the Hice University Instrument in 1968^ and 1974^, 

and HEAO A-4 in 1978^. No lines were seen at 511 keV in any of the 

nresolved line 
, respec- 
tively. These instruments employed scintillation detectors with 
energy resolutions at 511 keV in the 40-70 keV FVHM range. We 
report here the results of the first high— resolution gamma— ray 
observation of Cen A, made by an instrument employing Germanium 
detectors with a resolution of 2.2 keV FVHM at 511 keV. 


three observations, with 2o 
being 1.8 X 10“^. 8 x 10"^, 


flux upper limits for an i 
and 6.5 X 10” ph cm”^ s” 
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The GSFC Low Energy Gecasa-ray Spectroceter (LEGS) pecforcs 
hlgh>reoolutloo opectroscopy between ~ 70 keV and 8 MaV uelng an 
array of three cooled hlgh-purlty Gc dctectora surrounded by an 
active Hal scintillation shield. The detectors have an active voluae 
of 230 co’ and a peak effective area of 35.5 cd^ at 130 keV. 

The average ln~f light energy resolution rises froo 1.8 keV FWIM at 
70 keV to 3.5 keV at 2.6 KeV, At 511 keV the effective area Is 13.3 
cn^, and the resolution Is 2.2 keV FWHM. The active Hal shield Is 
~ 12 ca thick and colllcates the fleld-of-vlew of the detectors to 
~ 16* FWHM. The Instruaent la balloon-borne and Is oounted In a 
servo— controlled gondola that uses an altazlcuth pointing systea 
under nlcrocoscputer control; the pointing precision la ~ 0.5 . A 
detailed description of LEGS Is given by Paclesas et al . . 

The lostrus^nt was launched froa Alice Springs, Australia on 
1981 Kovc'^bcr 19 and observed Cen A for 3 hours at an average line— 
of-slght ataospherlc depth of 3.5 g ca . The observation was 
divided Into 20-alnute Intervals during which the telescope was 
alternately pointed at the source and away froa the source for 
background detcralnatlon. The source flux was calculated by 
subtracting the average background level froa each source Interval, 
correcting for detector efficiency and ataospherlc attenuation, and 
suaalng the resulting residual fluxes. 

Figure 1 shows the observed opectrua of Cen A, calculated as 
described above. The best-fitting power law of the fora A(E/100 
keV)“° has A - 1.05 x 10"^ ph s~^ keV^ and a - 1.59. The joint 

902 confidence error Halts (xf. +4.6; reference 6) for A and a 
are shown In the Inset. The data are consistent over the entire 
aeasured energy range with a power law of photon Index - 1.59, 
showing no evidence for a spectral break. The observed values of A 
and a are both Interaedlate In the range of values aeasured pre- 
viously for Cen A (reference 4, and references therein). 

The data were searched for features In the spectrua with the 
result that no statistically significant narrow or broad lines were 
seen. The 2o flux upper Holt for a narrow (< 3 keV FWKM) 511-keV 
positron annihilation line Is 9.9 x 10"* ph ca"'^ s"^, calculated 
froa the source flux In a 4-kcV wide bln' centered on 511 keV. This 
Halt Is a factor of ~ 1.9 larger than the upper Halts at nearby 
energies due to the Intense Inetruaental background line at 
511 keV^. 


DISCUSSION 

The lack of a 511-keV line In the Cen A epectrua Ic of 
particular Interest In light of the 511-keV ealsslon that has been 
observed from the central region of our own galaxy (see, e.g., 
reference 8) . In order to conpare the nucleus of Cen A with the 
galactic center source, we calculate here the ratio of the 511-keV 
line luminosity (or upper Halt to the luminosity) to the luminosity 
of photons with energies > 511 keV. This second quantity Is not as 
accurately deteralned for either source as lower-energy Integrals of 
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th*ff spectra but Is leorr relevaac for the conparlson; for instance » 
the positrons responsible for the galactic center 5ll-koV oiai|oion 
are likely produced by yY Interactions of photons with > 
energies^. Also, for the galactic center region there U 
considerable source confusion in the low-energy gaGsna-ray cfiaiurc*- 
oents made to date, whereas at energies greater than a few hundred 
keV the nuclear source doainatea^ • 

Since no flux measureeents exist for Con A b«t%feen 1 fteV and 1 
GeV. an extrapolation of our measured spectrum (Figure 1) was 
required to obtain the Cen A luminosity at energies > 511 keV* We 
performed the extrapolation using the spectral form for the 
differential flux 


dF 

dE 


0.159 

♦ (E/2000 


-2 -1 -1 
photons CB ‘ B keV 


O) 


In Annlogy to the observed npectrum for 3C 273 • • The spectral 

form In equation (1) tends toward the measured speetrun below I HsV 
and Is consistent with the SAS-2 flux up|H?r limits^ between 35 and 
200 HeV. Based on equation (l). the Cen A luminosity at enerRles 
> 511 keV Is 2.7 X 10^^ ergs s" . Tt\e LEGS 2o upcer_l Imlt for 
narrow-line 511-keV emission of 9.9 x 10"* ph cb"^ s’* rIvcs a line 
luminosity upper limit of 2 .4 x 10 erRS s . The upper Holt for 
tlie 1 Ine-to-cont InyuB ratio Is, therefore, 0.090; a ofollar analysis 
using the UEAO A-4“ line Unit and spectnen Rives a rntlo upper 
Holt of 0.073. 

For the galactic center, we use the HEAD C-l observations In 
the Fall of 1979. Depending on the assumption of the shape of the 
spectrum, the > 511 keV luminosity Is In the range 1.4—1 .7 x 10 
ergs s”', based on data extending above I MeV (G. Rlegler, private 
communication. 1983). The measured narrow-line 5ll-keV luminosity 
at that time was (1.8 t 0.2) x 10^' ergs s"* ‘ \ yielding a llne-to- 
contlnuura ratio In the approximate range 0.10-0.13. Both 
and continuum fluxes from the galactic center are variable 
the line Hux falling below detectable levels In recent taeasure- 
raents However, there Is evidence for a correlation betwen the 

line and continuum flux levels^**'" so that the I Inc-to-cont Inuum 
ratio of -0.12 may apply to other times than the Fall 1979 
measurement. Both the LEGS and HEAO A-4 2o upper limits for the Cen 
A lliie-lo-eoMtluuum ratio fall below the measured galactic center 
ratio. Their combined weight gives evidence that the emission 
regions In the nuclei of the two galaxies are different. This 
result Is not surprising given the factor of '• 10 difference 
in gamma-ray luminosity of the two sources. The higher temperature 
and activity level that one might expect In the nucleus of Cen A 
could produce a broadening of any 511-keV line emission. None of 
the observations to date have detected a broadened annlhl lat lon- 
radl.'itlon feature In the spectrum, hut the sensitivity tor line 
detection decreases as tlie line wldtli increases aK've tite 
Instrumental resolution width. 
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C«n A wae observed on 1931 Noveober 19 «Uh e ‘'^Ch-recolutlon 
gftoae-roy epectroacter. The photon opcctruo between 70 end SOg keV 
uan found to be well fit by a power lew of the fora 1.05 x 10 
(E/lOO keV)"^‘^’, with no evidence of a break. Ho llnoo or fenturee 

were oeen In the epeotrua; the 2a flux upper ® 

3 keV) 511-keV pooltron-annlhllatlon line la v.v x lu 
ph J^ l< A co^arUon with the observed 5ll-keV line fjuj fron 
the gnl.ctlc center In the Fell of 1979 tndlceten thet Cen * 

“nlL. thnn the flectlc tenter In 511-teV nnrron-lln. enl.elon 
relative to the > 511 keV contlnuun ealsalon. 
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PAIR PRODUCTION tlEAR THRESHOLD IW PULSAR MAGNETIC FIELDS 

A. K. Harding 

NASA/ Goddard Space Flight Center, Grecnbelt, HD 20771 

and 

J. K. Daugherty 

University of North Carolina, Asheville, NC 28804 

ABSTRACT 


In pulsar polar cap models, curvature radiation y-raya produce 
e'^e" pairs in the strong magnetic fields near the surface of the 
neutron star. While these Y“^nys have energies mc^, they also 

propagate at very small angles to the field, ouch^that the threshold 
condition, Ry> 2ncVoinO is Just barely satisfied when they pair 
produce* Threshold effects on the pair production att^enuation co- 
efficient, which are due to the discreteness of the c'^e* Landau 
states, must therefore be considered when computing the mean free 
paths of curvature radiation photons in pulsar magnetic fields* 

These effects, which are not Incorporated in the asymptotic 
eiqpression for the attenuation coefficient, have some interesting 
consequences for pulsar models, .^ince pair production is suppressed 
near threshold, the photon c^an free paths are lontrer than 
previously thought. In magnetic fields > 6x10^^ G, the pairs tend 
to be produced in the ground state Landau level and will not 
synchrotron radiate. Since synchrotron radiation is an essential 
ingredient in the electromagnetic cascades which produce low energy 
pairs above the acceleration region, pulsars with very high magnetic 
fields may not produce many pairs. 


INTRODUCTION 


The production of e'*‘e"* pairs plays a central role in most 
current pulsar models. Polar cap voltage drops which accelerate 
particles to ultrarelativiotic energies are limited by pair 
production discharges, ^ which continue in the form of electro- 
magnetic cascades above the acceleration region*^ These nair-photon 
cascades can generate large numbers of e^e** pairs which are thought 
to be essential for the coherent radio emission process. The most 
important mechanica for producing pairs near the polar cap is pair 
production by single photons in the intense magnetic field of the 
neutron star. Recent theoretical study of this process has provided 
a description of the behavior of the photon attenuation coefficient 
near threshold and the energy distribution of the pairs. In this 
paper, we discuss the implication of these results for pair 
production and electromagnetic cascades in pulsar magnetospheres. 


PHOTON MEAN FREE PATHS 

Magnetic pair production has been studied extensively, and 
almost exclusively, in the asymptotic limit of low fields, B«6 
and high photon energies » 2mc2, where B^^ = 4.414 x 10^ ^G is 
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the critical field otrencth. If ve define the unite «’5 b/2q and^B*^ 
= B/D Ch ■ c - 1), then this Holt takea the forn, C s . 

The pa?aLtcr C lo related to the nuabar of encrcctlcally 
Landau otatca available to the e o pair. In thin llnlt, the 
polerltatlon-averascd photon attenuation coefficient In the center 
of B 38 s“ fraoe, where the photon propacateo perpendicular to B, 

Bj^(«',B') ~ 0.23 I B* exp (^) evT^, X«l 

where y h u' .B'. In on arbitrary frane, where the photon propagatca 
at an angle ?^to B, Ite energy can be found by Lorents 
along B with velocity, B - coaB, and the reoult la « - u_j/eln0. 

Slnllorly, the attenuation coefficient becocea R - aln 0 BcM* 

To produce a pair, the photon center of oaaa energy iniot exceed 
threohold energy and the attenuation coefflcent cuat be non- 
negllglble. Theoe two conditions are Independent and can 
expressed In the following simple forms. 


(A) Uqj ■ w aln9 >_ 2 m 

(B) X > 0-1 

THRESHOLD 



where the latter cornea from 
the exponential dependence 
of R(oj', B'). Figure 1 
Illustrates where these 
conditions are oatlofled for 
different field atrengtho. 

In pulsar magnetospheres, 
high energy particles produce 
curvature radiation photons at 
very small angles to the field 
such that uslnO < 2n even 
though 01 2a. In order to 

pair produce, then, these 
photons taist propagate In a 
straight path until they 
acquire an angle to the curved 
field lines which satisfies 

(A) . In figure 1, this Is 
equivalent to moving upward 
along the diagonal lines of 
constant field strength. If 

(B) Is satisfied before (A), 

as it is for B > 0.1 - 4.4 

X 10^^ G, then the photon will 

pair produce very near 
threshold, where Eqn (1) lo no 
longer accurate* It has been 


Fig. 1. Pair production parameter x vs. "center of mass" photon 
energy for different field strength8(dlagonal lines). 
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found (Ref* S) that the quantum offocta of the dlecroto pair statee 
auppreas the pair production attenuation coefficient bolo^ the 
aayoptotlc Halt, auch that a reoro accurate e^rooslon la: 




exp (- 


-A f 




■35T 


-I 

CB , 


CM - 


> 1 




« I + .42 


CM 


-2*7 «,-.0038 

p 


( 2 ) 


The function ((o>\ B') approxlaetea the behavior of the exact 
attenuation coefficient near threshold, after averaging over the 
sawtooth pattern of cyclotron resonance spikes (cf> Ref« 6)* 
Since sln6 Increases with the photon path length, s, 
approxlnately as s/p, where p Is the rcdlue of curvature of the 
aagnetlc field, the isean free paths of the photons, froa either 
condition (A) or (B) and Eqn (2) will be 


X 

P 


O.l f (u* 


CM’ 


u)'B’ 



B’) 

, B’ < O.l 

B* > 0.1 


(3a) 

(3b) 


Due to the threshold condition, the isean free path Is constant above 
~ AxlO^^ G. Using Eqn (1), the approxlnate mean free path would be 



P 


0.1 

(i)*B“' 


(A) 


so that, 

^ « I + 210 B’ B’ < O.l (5a) 

*A 

10 B', B* > O.l (5b) 


The actual moan free paths of curvature radiation photons are 
therefore larger than those derived using the veil knovn asymptotic 
limit and this discrepancy Increases vlth field strength. 

The longer mean free paths will cause the voltage drop at the 
polar cap to be somewhat higher than previously catlmatcd^*^ 
although this should not be a large effect* The secondary particle 
yields of the cascades above the acceleration region will be much 
more sensitive to changes in the mean free path because the 
multiplicative effects of large numbers of photons are Involved* 
Longer mean free paths In a magnetic field which falls off with 
dtscance from Che star should decrease the number of pairs produced 
per primary particle In these cascades* If proposed effects due to 
the vacuum index of refraction in a magnetic field are present » \ 
could be even larger. ® 
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DISCRETE e‘*‘e~ PAIR STATES 

When the pair is produced^ the electron and positron nust 
occupy discrete energy otates with veil defined energies 
perpendicular to the oagnctic field* Tlie spacing of these states In 
energy increases with B, reaching AE -- oc^ at B . For a given 
photon energy and field strength, there Is a set of kineoatically 
allowed otates into which the pair can be produced* Figure 2 shows 
the number of these states as a function of photon energy and 
aagnetic field for curvature radiation photons which ate emitted 
parallel to B and travel one mean free path before producing the 
pair* At each magnetic field, the photon's energy must exceed a 
certain value before a pair can be produced with at least one member 
in an excited state. In fields > 5 x 10^^ G, a significant fraction 
of photons in the curvature radiation spectrum (for typical pulsar 
parameters) will produce pairs in the ground state* 

Models of pulsar cascades^ have shown that synchrotron 
radiation Y*rays from secondary particles are necessary to sustain a 
cascade with several photon generations and high pair yields* 
Curvature radiation from the pairs Is much less efficient and unless 

Ep 


2xt0’*ev 5xt0'*ev 10’*«v 



Fig* 2* Number of energetically allowed e'^e* states vs. curvature 
radiation photon energy for different magnetic field strengths with 
constant radius of curvature, oy = p/lO^ cm. The top scale is the 
primary particle energy E for which the critical energy * 

(3/2) (Ep/mc2)3 c/p - 0 ). 

radii of curvature are significantly less than dipole, the cascade 
terminates after one photon generation. In magnetic fields high 
enough to suppress production of pairs in excited states, 
synchrotron radiation cannot take place (the photon energy is fed 
into particle motion parallel to the field). The efficiency of the 
cascades should therefore decrease rapidly above some field 
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strengtht depending on the crieleel frequency a * (3/2) 

(Ep/nc^)^ c/p of the prlnaty particle curvaturc^opectrua (Figure 
3). Cascade efficiency (ie. the ratio of secondary to prinary 
particles) will increase with field strength below this lioitt 
however. For values of the primary particle energy E predicted by 
various acceleration codclo copious pair production by 

cascades near the polar cap would not be expected to occur for field 
strengths above B ^ 6-8 x 10 ^^ 



CONCLUSIONS 

This 

reevaluation of pair 
production by 
curvature radiation 
photons in pulsar 
magnetic fields has 
shown that threshold 
effects on the 
attenuation 
coefficient can have 
significant 
consequences for 
pair-photon 
cascades. In 
particular, there 
should be a maximum 
in the pair yield per 
primary particle as a 


Fig. 3. Primary particle energy (l.e. polar cap acceleration 
voltage) vs. magnetic field strength above which all curvature 
radiation photons (with m j< 3 « ) produce pairs in the ground 

state. 


function of magnetic field strength with the yields diminishing 
rapidly in fields > 6 x 10 ^ inhere synchrotron radiation is 
suppressed. Since the copious production of e"^e" pairs is necessary 
for coherent radio emission in most pulsar models, this may imply 
that there is a high magnetic field cutoff for radio pulsars. 
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AIRIIHILATION AND Ttf£ COSmC X-RAY BACKGROUND 

Demosthenes Kazar.^s cad Rtchard Arrlck ^afer 
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and 
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University of Maryland 

ABSTRACT 

The possibility that the processes responsible for the Cosmic 
X-ray Background (CXB) would also produce an e^-e"^ annihilation 
feature is examined* Under the assumption that these processes are 
thermal, the absence of a strong e'"-e^ annihilation feature places 
constraints on the compactness (L/R ratio) of these sources* 
Observations favor sources of small compactness ratio* 

INTRODUCTION 

The fact that the X-ray sky is dominated by an isotropic 
component (the so called Cosmic X-ray Background, hereafter CXB) had 

been established by the earliest 
X-ray astronomy observations^. 
The subsequent satellite X-ray 
observations, especially by the 
A-2 and A-4 experiments on HEAO 
allowed the detailed 
spectral determination of CXB. 
The observed spectrum in rhe 
region 3-150 keV, along with the 
higher energy data is shown In 
fig 1. The HEAO 1 experimenters 
have found that thermal brems- 
strahlung from an optically thin 
plasma of temperature 40 ± 5 keV 
provides a remarkably good fit 
to the data from 3 to 100 keV. 
Interestingly enough, no studied 
» to' lO* 10 * to* 10 * io‘ populat ton of sources is known 
e«ERCv (uvi have a thermal spectrum with 

the required properties - 

Fig. I - The unresolved X and garoira-ray 
background. (From refs 2»3*12*13) 

One can of course contrive to combine sources with a variety of 
spectra emitting over a range of red shifts to produce the observed 
total background spectrum, even If the Individual spectra arc 
different than that of CXB^* . The shape of of the spectrum however 
clearly suggests a thermal dlstrlbut l*. i of ratlu^r specific 
temperature and It would be more natural If the CXB could be 
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•explained ae such« 

A thin thenaal breesBtrahlung from a heated Intergalactlc 
medium will provide such a spectrum However* the energy 

ret:juired to heat a diffuse uniform taedium tc euch a temperature is 
quite large. Also* the low densities end consequent long electron- 
ion coupling times provide difficulties in attaining and maintaining 
a Maxwellian distribution for the medium^. Clumping of the medium 
would reduce the input energy requirements and the coupling time 
scales. Taken to extremes* the clumps might be reduced to a size 
comparable to galaxies* or smaller* becoming •'compact'* sources of X- 
rays. 

For cither the heated intergalactic medium or the "compact 
source" models for the CXB* the bulk of the emission originates at 
redshifts -- 2-3 (or even larger) so that the corresponding source 
temperature would be kT > 100-200 keV. For these temperatures a 
sufficiently compact source will produce electron-positron pairs 
from the tails of the photon and particle thermal distributions. 
Under certain conditions the positron abundance would be sufficient 
to produce an observable e*^— e annihilation feature in the CXB. 

THE POSITRON ABUNUAl^CE 

In a thermal plasma of temperature kT > 100 keV it is possible 
to produce positrons at significant abundances by ee* ey or yt 
collisions since a non -negligible fraction of the particles and 
photons at the tails of the distributions fulfills the pair 
production threshold condition. Their steady state abundance Is 
determined by the balance between pair production and annihilation 
reactions (for a detailed treatment see ref. 8). 

In the cases of interest for the CXB the dominant positron 
production is due to yy reactions so the ee, reactions will not 
be considered further . (Their cross sections are correspondlrgly 
smaller by and a). The approximate expressions for tht production 
and annihilation the rates are: 


V * hv V V • 


( 1 ) 


8 ^+'’- 


where n • n * are the positron and electron number <*enslties, 
the T^iomson cross section and x * 2 m c^/kT. f(x) is a function of 
the temperature only, resulting from the averaging of the photon- 
photon pair production reaction rate over the photon distribution 
functions- It is given by 


f(x) - 


-X 


^ 

x(x-K:^)* 


where * 1.1A3 and C 2 * 3.63 are constants.KY Is the normalization 
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of lUit breo^strohlani;: photon disttltMit ion function detoroincd by tho 
condition 


• % -E/kT 
o 


E dE • or 


•^v . kUt!> . 

^ kT 

where L* R, t* T ore respectively the luoinosity* r«dius« optical 
depth and teoperature of the eourcet. The steady state condition 
Rives 

tt ^ 


5,- '<*> - 


(n^^n_>* 


Since the cowpset sources presently considered are prcsuasbly 
gravi tat ionsl ly bound, one can scale L and R by their |:ravi t at iona 1 
units i .e • the Eddlnitton luoinoslty (L^ • 10^^ M erj; s^^) and the 
Schwarrschlld radius (K^ - 3 10^ M c») where H is the oass ol the 
sources in solar niasses. For this we Introduce two parameters F and 
f, both with cKpected values between 0 and I, defined as 

F I U\^ and f 5 R /R . 
fc s 

Their product Ft * L/R is a mass Independent measure of their 
compactness, A hi|:h value of FI corresponds to a hiRh density of 
photons and hence e^e pairs. Hie additional factc>r ol R in R can 
be espressed In terms ot t»ie optical depth t, since {t\ ) o^ R • 

% St where s is the aspect rat lo of the source, t .e . the ratio of 
its lary.est to its shortest dimension. We can therefore reespress 
eq (3) in terms of the positron abundance X i n^/n by: 


2x\0 FJ 


%/2 X 

> 


1/? r»l ^ 

*st * 14V 


roMPAKisiW n> onsKKVArioN 

hijua: ion (i) can be dltecllv related lo observal Ions ot the 
CXb • l^t A be t tte rat to %>t t tie antiihl lat ion to I lie bt emst t atiluitu; 
spectral luratnosiiv at t tie peak enei ot the anni tit 1 at ton leature. 
K • m 4 kT. Since no obvious annthtlatton teatuie is observed in 
t lie A V 1. I’siny, t lie tesults o! ret. 'J. one can relate 

V to A K t'V 
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Eq (6) is Independent of Che 
caass of the sources and Holts 
(as a function of the redshift Z 
at which the CXB was produced) 
their L/R ratio to be coapatiblc 
with the absence of a proninent 
annihilation feature in the 
CXB. The optical depth t of the 
sources is unknown, however it 
is constrainted to be t £ 3 
otherwise the corresponding 
self-Conptonlzatlon Wien peak 
should be apparent in the 
spectrua^®. Another unknown 
paraceter of Che sources is the 
aspect ratio. For spherical 
sources s* I* although it could 
be s»l for sources of thin disk 
geocetry. It appears, however, 
s » 1 that for sources ealting 
close to their Eddington 
lunlnos t ty ^ ^ • The constraints 
ieposed by eq (6) are shown in 
fig 2. For exaoplc, the absence 


Fig. 2. Linits on source lualnoslty F 
(in units of Lg) vbs redshift of CXB 
production, z, for different source effective 
sizes, s/f (in units of R^.). The tcaperature 
is the source teoperature in the ennitter frame. 


of a prominent annihilation feature at 100-150 keV, constraints 
compact sources (.i^O.l) operating at to emit at a small fraction 

of their Eddington lualnoslty. Due to the absence of positive 
detection of such feature all the derived constraints are in fona of 
Inequalities which however nay prove Important In understanding the 
nature of the sources of CXB. Fig 3 shows the superposition of the 
energy spectrum of thermal bremsstrahlung with T - 200 keV and an 
annihilation feature of the same temperature and A ■ 1. Such a 
feature appears detectable, though unfortunately In this example It 
peaks in the observer frame energy range E > 100 keV, which Is 
probably dominated not by the thermal component but by the 
contribution of active galaxies. This Is currently very poorly 
known and we may always have to be content with the existing upper 
llmts. A positive detection an annihilation feature In the CXB 
spectrum, rather than the above qualitative upper 

provide an Important contribution to the understanding of the CXB s 
origin. Such a detection will determine precisely the emitter-frame 
temperature (and hence the redshift) of the sources. It will signify 
that they are compact and It will provide a consistency check for 
the temperature of the ternal component. If thU feafir.j Is 
detectable in the CXB spectrum, as a whole It may someday be studied 
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Ln individual spectra of compact 
objects • if ouch objects cake up 
the background* However ouch 
studies will probably require 
high energy imaging capabilities 
beyond those planned for the 
next generation of expcricents. 


Fig. 3. Superposition of brensstrahlung and 
annihilation radiation energy spectrum for 
A • I and T • 200 keV. 
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ABSTRACT 

Observations of the e*-e" annihilation radiation from the Galatic Center suggest that 
something truly extraordinary is occurring there. We review the observations of this 
intense, time-varying, 0.511 MeV emission and discuss the implications of these and other 
recent observations on the positron production process, the annihilation region and the fun- 
damental nature of the Galactic Center source. 


INTRODUCTION 

A year ago we reviewed* the possible origins of the e+-c" annihilation radiation from 
the Galactic Center and concluded that the most likely process for producing the annihilat- 
ing positrons is photon-photon pair production in the vicinity of a massive black hole at the 
center of our Galaxy. We briefly summarize here the evidence and arguments leading to 
this conclusion, and at the same time discuss the implications of new observations and cal- 
culations that shed further light on the problem. 

We briefly summarize th® observations^ and then discuss in turn their implications on 
the nature of the annihilation region, the positron production process and the Galactic 
Center source itself. 


OBSERVATIONS 

Intense positron annihilation radiation at 0.511 MeV has been observed from the 
direction of the Galactic Center for over a decade. This emission was first seen in a series 
of balloon observations^'* with low-resolution Nal detectors starting in 1970. But it was 
not until 1977 that the annihilation line energy of 0.511 MeV was clearly identified with 
high-resolution Ge detectors flown by Leventhal, MacCallum and Slang*. The latter obser- 
vation also revealed that the line is very narrow (FWHM ^ 3.2 keV) and suggested that 
the continuum telow 0.511 MeV may include a significant contribution from three-photon 
positronium annihilation, consistent with ~ 90% of the positron annihilation taking place 
through positronium formation. 

The existence of this very narrow line was confirmed by Riegler ct al.’ with Ge detec- 
tors on KEAO-3 in the fall of 1979. These observations set an even more stringent limit 
on the line width (FWHM < 2.5 KeV) and determined the line center energy as 510 90 ± 
0.25 keV. 

The HEAO-3 observations^ also provided new information on the location and spacial 
extent of the emission region and most important showed that the line intensity varies 
significantly in time. In particular, these observations showed that the line emitting region 
is smaller than the angular resolution of the detector (35* FWHM) and that the direction of 
the source coincides with that of the Galactic Center, within the observational uncertainty 
of ± 4“. Moreover the observations showed that the 0.511 MeV line intensity decreased 
by a factor of three in six months, from (1.85 ± .21) x IQ--* photons/cm^ sec in the fall of 
1979 to (0.65 ± .27) x lO"* pholons/cm^ sec in the spring of 1980. This variability has 
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been conflnncd by b:lloon>bome Ge detector ob^rvetions**’. Observations'^ vrith a Nal 
detector in the fdl of 1977 could elso indicate a variation on a time scale cs short as 
10 days, but it seenu much more likely that the hi^er 0.S11 MeV intensity observed with 
this detector results from a terser diffuse galactic component seen in its much greater 
(KW'FWHMieeldofview. 

Observations of continuum emission in the hard X*ray and gamma*ray bands have 
recently been reviewed by Matteson". The hard X-ray emission is also time variable and is 
weakly correlated with the variability of the O.Sll McV line '(e.g. Ref. 9). These observa- 
tions set an upper bound of 2 x 10^^ erg/sec on the Galactic Center continuum luminosity 
at photon energies > nicC^, since only part of this emission may come from the same 
source as the annihilation radiation. 

The luminosity of the Galactic Center region at various photon energies, implied by 
these and other observations are summarized in Figure I from Reference 1. 



Rgure 1. The luminosity per unit InE as a function of photon energy. E. from the region 
around the Galactic Center. Data are shown for the compact (< 10'^ cm) nonthermal ra- 
dio sourcc'^ the — 3pc dust ring'^, the nonthermal infrared source IRS 16 (Ref. 14), the 
soft X-ray emission (< 3pc) from the EINSTEIN satellite measurements'*, the hard X-ray 
emission from HEAO-I (Ref. ID, the 511 KeV line and positronium continuum*’-’, and the 
gamma ray emission from HEAO-1 (Ref. 1 1) and COS-B (Ref. 16). Also shown as dashed 
curves arc the blackbody luminosity at ~ 31,000 K required" to account for ionization in 
the warm IR c'ouds within < Ipc of the Galactic Center and a blackbody luminosity at 
120 K as inferred" from the far infrared observations of the •— 3pc dust ring. 
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THE ANNIHILATION REGION 

The nature of the positron annihilation region is constrained by the intensity varia* 
tions, the line width and the line center energy. Tlie size of the region should not exceed 
about 10“ cm, the distance traveled by relativistic positrons in 1/2 year. The density of the 
gas in which the positrons annihilate should be larger than 10^ H/cm^ the minimum den- 
sity required to slow them down in 1/2 year, but less than 10“ H/cm* in order not to break 
up the triplet positronium before it annihilates. Since triplet positronium could also be bro- 
ken up (R. McCray, private communication) by photons of energy > 6.8 eV, the energy 
density of such photons should not exceed ~ 10^ crg/cm* in the annihilation region, or if 
the 31,000 K emission” (Rg. 1) comes from a single source the annihilation region must 
be > 2 X 10 “ cm away from it, assuming a positronium ionization cross section of 
3 X 10~“ cm*. 

The observed line width requires** that this gas also be at least partially ionized 
(n, < O.ln). If the gas were neutral, the line width would be larger than observed because 
it would be Doppler broadened, not by the thermal motion of the gas, but by the velocity 
of energetic positrons forming positronium in flight by charge exchange with neutral hydro- 
gen. In a partially ionized gas, however, positrons lose energy to the plasma fast enough 
that the positrons thermalize before they annihilate or form positronium. The line width 
thus reflects the temperature of the medium, so that the observations require a temperature 
5 X 10^ K. The line width further limits any velocities of rotation, expansion or random 
motion to < 700 km/sec, while the line center energy implies a bulk velocity along with 
line of sight -90 < v < +200 km/sec and a gravitational redshift z < 7 x 10^ 

The strongest of these constraints are summarized in Table I. 

Table I 


CONSTRAINTS ON THE e+-e" ANNIHILATION REGION AT THE GALACTIC CENTER 
Physical Parameter Constraint Observation 


Size 

Gas density 
Ionization state 
Temperature 
Rotation, expansion 
or random motion 
Bulk motion along 
line of sight 
Gravitational redshift 


< 10'* cm 
> 10* H/cm* 

tte/n > 0.1 
£ 5 X I0< K 
< 700 km/sec 

-90 < V < +200 km/sec 

z < 7 X 10< 


variability 
variability 
line width 
line width 
line width 

line center energy 

line center energy 


As we previously suggested” possible annihilation sites which satisfy these constraints 
are the warm clouds** and the compact source IRS 16 (Ref. 12), observed within the cen- 
tral parsec of the Galaxy. 


THE POSITRON SOURCE 

The nature of the positron source is also strongly constrained by the observed varia- 
tion of the 0.511 MeV intensity and by observations at other wavelengths. The decrease of 
a factor of three in the line intensity in six months clearly excludes any of the multiple, 
extended sources, such as cosmic rays, pulsars*®, supernovae*', or primordial black holes , 
previous proposed. Instead, it essentially requires a single, compact (< 10'* cm) source 
which is apparently located cither at, or close to, the Galactic Center and which is 
inherently variable on time scales of six months or less. 
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The observed 0.51 1 MeV line intensity of — 2 k 10“^ photons/cm* sec requires at the 
distance of the Galactic Center (— 10 kpc) a positron annihilation rate of 4 x 10*^ tVsec, 
if ~ 90%^ of the positrons annihilate via positronium. This rate corresponds to minimum 
luminosity of ~ 6 x 10^^ erg/sec in both the line and thrcs>photon continuum. With such 
a luminosity the Galcctic Center is the most luminous gemma>ray source in the galaxy. 
The uniqueness of this source makes it unlikely that it results from the chance occurrence 
of the youngest supernova or pulsar along with line of sight to the center of the galaxy. 

Tlie strongest constraints on the various positron production processes are set' by 
observations of continuum emission at energies > m«c^ from the direction of the Galactic 
Center". When compared with the annihilation radiation luminosity, the continuum 
gamma ray luminosity implies a very efneient positron production process, one in which 
more than 30% of the total radiated energy > goes into electron-positron pairs. If the 
positron production occurs on time scales comparable to that of the observed variation and 
in an essentially optical thin region which emits isotropically, only photon-photon pair pro- 
duction can provide the required high eflkiency. 

We considered* two geometries for the positron production region: a spherical 
volume in which e'*'— e~ pairs are produced by photons interacting isotropically and a beam 
in v/hich the pairs are prc^uced by photon interactions only at small angles. 

The most efneient pair production occurs in isotropic interactions of photons at ener- 
gies close to nicC^. The pair production rate Q in a spherical source of radius r may be 
approximated by 

Q ~ !/4 <orc> r^ 


where <orc> is the average pair production cross section times the velocity of light, equal 
to ~ 3 X lO"'* cmVsec for black body photons of temperature ~ mjC^ end n^ is the pho- 
ton number density. Assuming that the source is optically thin, the photon density can also 
be related to the continuum luminosity at energies mec' by 

«nyC4irr' 

L~ j . 


where < is the average photon energy and r/c is the photon residence time. Combining 
these two equations and setting e — meC*. we see that for a given continuum luminosity 
the positron production rate depends only on the source size, such that the radius. 


r 


3 <gc> 

8wc* (mec')' 


|2 

^-6x10 


r« (cm). 


From the observed luminosity limit of L 2 I x 10^* erg/scc and a production rate Q 
equal to the annihilation rate of 4 x 10^^ eVscc, the radius of the positron source must be 
6 X 10* cm. 

Pair production by isotropic photon-photon interactions thus requires an exceedingly 
compact course, but with a high luminosity. The most obvious candidate is a blackhole. 
But if this source is a blackhole releasing gravitational energy of accreting matter close to its 
Schwarzschild radius, then it must have a mass < 10^ M„ which is much smaller than the 
masses of 10^ to 10’ M* blackholes that have been suggested*’’’’ at the Galactic Center. 
Yet such a small size would be consistent with arguments by Ozernoy’^ that the Galactic 
Center cannot contain a blackhole larger than about 10’ M„ if tidal disruption of stars is 
the principal source of the accreting matte, on which it grows. 

The photons needed to produce the pairs could themselves be produced in a hot 
accretion disk around the blackhole”. A luminosity of 2 x 10’® erg/sec requires an 
accretion rale of -- 3 x 10"® M^yr which would form a — 300 hole in the age of the 
Galaxy. A major fraction of the e"*" pairs produced by photon-photon collisions above the 
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disk could then escape from the source region before they annihilate, a constraint (Tsbte 1) 
set by the absence of any measurable redshifl in the energy of the ennihiletion line. 

We turn now to the cltcmative geometry of pair production by small cngle photon 
interactions in a beam, which may be produced*^” by dynamo action in a msgnetic field 
accreting onto a blackhole. We previously showed' that with a beam the constraint on the 
size of the production region could be greatly relaxed, but at the expense of a much higher 
beam luminosity in gamma rays of energy > > m«c^ This possibility for producing posi* 
trons in the Galactic Center was first suggested by M. Bums (private communication 1932) 
and differs from Novikov’s^^ model (in this volume), which relies on the relatively less 
efficient production by beam photons interacting with gas in a cloud. 

The pair production rate Q for small angle (0 — rt// ) photon interactions in a beam 
of radius rt, and length /may be approximated by 

Q — '/4 n,* <<rVi> irr^V, 

where v^ — (rt//)c is the mean transverse velocity of the interacting photons, end n^ is the 
density of those photons with energies greater than the small angle pair production thres* 
hold (f ‘fhis density can be related to the beam luminosity of such pho> 

tons by 

Lt — U/r^ m«c* n,c irrj,*. 


Combining these two equations, we see that the beam radius is 


fb- 


<<rc> 


2irc* (m.c»)* Q / 


(-v)* - 8 X l0-« ^ (•v)*cm. 


Thus for a pair production rate Q of 4 x 10'^ eVsec the beam radius could be as big 
as — 10" cm, or equal to the Schwarzschild radius of a 3 x 10* blackhole, if the beam 
luminosity at photon energies greater than 2S MeV were os high as half the Galactic Center 
bolometric luminosity limit of 3.5 x 10*' erg/sec, and the aspect ratio of the beam were 
0.02, corresponding to angle of 1*. 

The resulting pairs would also have energies of ■>- 25 MeV, comparable to those of 
the photons which produced them. But they could be stopped end annihilate to give narrow 
0.511 MeV line emission, if the beam hit a gas cloud. The bulk of the pair energy, 
amounting to — 10*** erg/sec, would be dissipated in heating the gas which could in turn 
rerodiate it isotropically os thermal radiation consistent with the constraints on the 
30,0(K) K luminosity. Since the radiation yield of ~ 25 MeV electrons and positrons is 
3%, their bremsstrahlung could also be consistent with the hard X«ray and gamms'ray 
luminosity limit of 2 x 10" erg/scc. 

The detailed energetics of both of these geometries, however, are still under study. 


SUMMARY 

The observed time variations and line width of the e*-e" annihilation radiation from 
the Galactic Center require that the positrons be produced essentially by a single source and 
that they annihilate in an ambient gas of density > 10' H/cm', ionization fraction > 10%, 
temperature < 5 x 10* K, and confined to a region of size < 10'* cm. Such conditions 
may exist in the warm clouds and the compact source IRS 16 within the central parsec of 
the galaxy. 

The limits on the accompanying continuum emission at energies > m«c* set strong 
constraints on the positron production process, requiring an exceedingly high efficiency, 
such that > 30% of the total radiated energy > m«c' goes into e*-c~ pairs. The most 
likely mechanism appears to be pair production in photon*photon collisions in the close 
vicinity of a massive blackhole. either near the hot (kT — m,c*) inner j«rt of an accretion 
disk around a — 10' M, blackhole, or in a beam of — 25 MeV photons produced by a 



ORIGINAL PAGE IS 
OF POOR QUALITY 


— ICi* M, hole. In cither case the absence of any measurable redshiA in the line center 
encrijy requires that a Icrne fraction of the positrons escape from the central source end 
annihilate at great distances from the hole (> 10^ times the Schwarzschild radius). 

ACKNOWLEDGEMENTS 

The work of R.EL. was supported by NASA grant NSG-7S41. 

REFERENCES 

1. R.E Lingenfelter and R. Ramaty, The Galactic Center (Am. Inst. Physics, New York, 
1982), p. 148. 

2. C.J. MacCallum and M. Leventhal, these proceedings (1983). 

3. W.N. Johnson, F.R. Hamden and R.C. Haymes, Ap. J., 172 LI (1972). 

4. W.N. Johnson and R.C. Haymes, Ap. J., 184, 103 (1973). 

5. R.C. Haymes ct al., Ap. J., 201, 593 (1975). 

6. M. Leventhal, C.J. MacGallum and P.D. Stang, Ap. J., 225, Lll (1978). 

7. G.R. Riegler et al., Ap. J.. 248, L13 (1981). 

8. M. Leventhal ct al., Ap. J., 260, LI (1982). 

9. W.S. Paciesas et al., Ap. J. 260, L7 (1982). 

10. D.M. Gardner et al.. The Galactic Center (Am. Inst. Physics, New York, 1982), 
p. 144. 

11. J.L. Matteson, The Galactic Center (Am. Inst. Physics, New York, 1982), p. 109. 

12. K.Y. Lo et al., Ap. J., 249, 504 (1981). 

13. I. Gatley, The Galactic Center (Am. Inst. Physics, New York, 1982), p. 25. 

14. E.E. Becklin et al., Ap. J., 219, 121 (1978). 

15. M.G. Watson ct al., Ap. J., 250, 142 (1981). 

16. 8.N. Swanenburg et al., Ap. J., 243, L69 (1981). 

17. J.H. Lacy ct al., Ap. J., 241, 132 (1980). 

18. R.W. Bussard, R. Ramaty and R.J. Drachman, Ap. J., 228, 928 (1979). 

19. R. Ramaty and R.E. Lingenfelter, Phil. Trans. R. Soc. Lond., A301, 671 (1981). 

20. P.A. Sturrock and K.B. Baker, Ap. J., 234, 612 (1979). 

21. R. Ramaty and R.E. Lingenfelter, Nature, 278, 127 (1979). 

22. P.N. Okeke and M.J. Rees, A. Ap., 81, 263 (1980). 

23. D. Lynden-Bell and M.J. Rees, M.N.R.A.S., 152, 461 (1971). 

24. L.M. Ozemoy, Large Scale Characteristics of the Galaxy (Reidel, Dordrecht, 1979), 
p. 395. 

25. D.M. Eardlcy et al.. Ap. J.. 224, 53 (1978). 

26. R.D. Blandford, Active Galactic Nuclei (Cambridge Univ. Press. London. 1979), 
p. 241. 

27. R.V.E. Lovelace, J. McAuslan and M. Burns, Particle Acceleration Mechanisms in 
Astrophysics (Am. Inst. Physics, New York, 1979), p. 399. 

28. N.S. Kardashev, l.D. Novikov, A.G. Polnarev and B.E. Stern, these proceedings 
(1983). 


62 





N83 27937 ^ 


CRITERIA FOR GRASAR ACTION IN ASTROPHYSICAL SOURCES 
J. M. McKlnlejr* and R. Ranaty 

ORIGINAL PAGE IS Laboratory for High Energy Astrophysics 

OF POOR QUALITY NASA/Goddard Space Plight Center 

Greenbeltt Maryland 20771 

ABSTRACT 

The theory of gamaa-ray aapllflcatlon through stlnulated 
annihilation radiation (graoar) was developed by P^uaaty, McKinley 
and Jones ^ (hereafter RHJ). For gccsia-rcy bursts slollar to the 
March 5, 1979 burst, an observed annihilation line of width <0,03 
MeV would Inply a grasar source* The alnla^ pair density needed 
for the onset of graoar action Is ^ 10^0 cm ^ and the peak of the 
grasar line, without a gravitational redohlft, la at < 0.5 MeV. 

INTRODUCTION 

An ecilsslon line at MeV, generally believed to be 

spontaneous, optically thin and gravitationally redshlfted e"^-c" 
annihilation, has been observed^* ^ from several gamma-ray bursts. 

The t^asured photon fluxes and the likely distances and sizes of the 
burst sources suggest^, however, that for at least some bursts the 
source regions are optically thick. Compton scattering and Y“Y pair 
production are the principal mechanises that would remove photons 
from an emission line In a gamma-ray burst source. 

In a detailed calculation of the emlsslvltlcs and the absorp- 
tion coefficients for two-photon pair production and annihilation 
and the accompanying Compton and Inverse Compton scattering, RIU^ 
showed that an emission line at '^J.43 MeV could be produced In an 
optically thick source without a gravitational redohlft by amplified 
annihilation radiation. In the present paper we consider the 
observational signatures thaj: would Imply the existence of such a 
grasar source, the minimum e’^-e** pair density needed for the onset 
of grasar action and the reasons that the central energy of the two- 
photon grasar line is at < 0.5 HeV. 

OBSERVATIONAL REQUIREMENTS FOR A GRASAR SOURCE 

The need for a graoar source can be seen from a relationship 
between the width of an observed annihilation line, AE, the fluence 
In the line, F, the duration of the emission of the line photons, 

At, the source distance, d, and the source projected area, A. We 
derive this relationship by comparing the minimum pair density 
required to produce the line by nonampllfied radiation (l.e. 
radiation with total absorption coefficient Kj>0) with the maximum 
pair density allowed by broadening due to pair degeneracy^ and the 
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oboecved upper Ilsit on the line width* We first calculate the 
ainimua pair dsnsity* 

The equation of radiative trenafer, dl/dt ■ J -R-1, requires 
that for nonanplificd annihilation radiation, the ratio of the 
spontaneoua annihilation eaissivity to the obaorptlon coeffici'^nt, 
j/K_, be larger than the radiation intensity I* Othenrioe 
dl/ 9 t <0 and no observable radiation Is prodveed* BHJ^ carried out 
detailed calculation of J and for pair placaas in a bath of 
aBbleot photons* Here ue estioatc the annihilation eaissivity, 

J, by a sinple expression which is consistent with their 
calculations, 

j « 2wr^c n^/(4w&B), (1) 

where n^. is the pair density (we assume equal positron and electron 
densities) and r “2.82x10“^^ cn. To calculate the olnlaua pair 
density, we can ignore the ambient photons* Thus > K^, where 

K- « 2n. vr^ (2) 

c + o 

is the Compton absorption coefficient* The intensity is given by 

I - . (3) 

^ AAtdE ' ’ 


By combining Eqo* (1)» (2) and (3) and using the condition for 
observable nonamplified radiation (J/K.j«>I) ve obtain 


"+ > Sr " ‘TS>^<r 




We next consider the upper limit on n^. For nonamplified 
annihilation, degeneracy broadening* sets a lower limit on the width 
of the annihilation line, 

AE >0*8. pc-0*8(3n2)^^^ he (6) 

where Pp is the Fermi momentum of the pairs* An observed line width 
or upper limit on the width, therefore, sets an upper limit on 

< 8*5x1027 ca“3 (AE/0*lMeV)3* (7) 

By combining Eqs. (5) and (7) we obtain a necessary condition 
for the production of an observed line by nonamplified annihilation, 

(ro) (ilk) > 

The most Intense astrophystcal annihilation line observed so 
far Is that secn^ In the spectrum of the 1979, March 5 burst whose 
source direction coincides^ with that of the supernova remnant N49 
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lo the Large !!agcllaolc ttoad CUiC). For this baret^ P - 7 j&otoas 
ca~2, AS <0.13 toV, and If the Itne-fonsatloa ties eqcala the 
duration of the Icpulslve phnse. At " 0.15 oec. Purthecrore, If ti^ 
eouree of the burnt la Indeed In the UCC, d ■ 55 l^pc* Ihcn frea SJq. 
(8), oonacpllflcd annihilation requiree that A > 3 ka * This 
condition be aatlofied on a neutron etar eurface elnce ita 
projected area io ~ 300 ka^. 

For the March 5, 1979 burnt* howercr, both At and AE cmld hare 
been algnlflcantly aaaller than the valuea given above. In 
particular* the line foraation tine could be Icon than the total 
Itapulaive phase duration* oince the annihilation tlc>e of the pairs 
is ertrenely short® (~ 10"IS sec). If future observations should 
indicate short duratlona and narrow lines* a strong case would eaist 
for grasar sources in astrophysics. 

TSEESBOLD DEHSlTf FOR GRASAR ACTIOO 


Spontaneous pair annihilation into two photons and the inverse 
process of two“photon pair production are necessarily accocpanlcd by 
stliBulated annihlUtlon (e.g. Bef. 1). The absorption coefficient 
due to stltaulated annihilation* Kg^* is always negative* but the 
generation of asplifled annihilation radiation (i.e. grasar action) 
requires that the total abcorptlon coefficient be negative for at 
least sose photon energies* 


0>Kj-Cc+Kpp+KsA • 


Here is contributed by pair production. is related^ to the 

spoot^eouo ealsolvlt^ of antilhllAtlon photons, j, 

( 10 ) 

where E is the photon energy. This expression* like its analogue 
for any other radiative process* can also be found directly frea the 
Einstein A and B coefficients. 



RMJ showed^ that for grasar action to occur the pair density n^ 
oust exceed a threshold value. To provide an estlaate of this 
threshold, we consider the case of a cold pair plassa (kT « Feral 
energy) with no aablent photons present. Such a degenerate systea 
has the lowest threshold density* since ^ I® reduced by 

the degeneracy. In this case we can use S^s. (1), l2), (6) and (10) 
and E * ac^ to obtain 



,fi .2 2/3 c . ,0 
JT7I f » 5.9x10 


-21 


Icn 




( 11 ) 


where f is the fraction of the degenerate positrons and electrons 
which can contribute to Coapton scattering at the photon energy of 
Interest. The threshold for grasar action occurs when ■*SA • 
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i«e* at a pair tenalty 

( 12 ) 

Accurate 

ealculetiea of and 
Kg for a desenerato 
posltroa-oLectroa placsa 
have beea carried out by 
IUU^« In Pig. 1 ve ohc9 
their resulto for a pair 
density just above 
threshold. As can be 
seen, belou the Feml 
energy Kq Indeed 
decreases rapidly with 
decreasing photon 
•jnergy. By cosparing 
the threshold density 
of ~103® cn obtained 
froa the nuscrlcal 
calculations, vlth tttat 
given In Eq. (12), ve 
see that f « 0.6, l.e. 
Figure 1. Absorption coefficients vs. photon 
energy In a degenerate positron-electron plasaa. 

approxlcately AOZ of Che degenerate positrons and electrons cannot 
contribute to Cocpton scattering at ~0.5 KsV. 

The threshold density for graoar action vould be lover than 

1030 

ca"^ If the Coc^tOQ scattering cross section were less 
than Indeed, in the presence of a strong oagnctic field, at 

photon rrequencies belov the cyclotron frequency, v , the Cos^ton 
cross section is substantially reduced* But the observation^ of 
cyclotron absorptloo features at 40*60 keV in cany gcasa-ray bursts 
euggests that for e"^-e annihilation photons E > 10 hv • Daugherty 
and Ventura^ have shovn that at such photon energies, 8he scattering 
cross section for photons directed parallel to the field is very 
nearly Altogether, it seems that the density of pairs in a 

grasar source taist erceed 10^® cn~^. 

THE CENTRAL ENERGY OF A GRASAR LWE 

Froa their numerical calculations RMJ^ found that the naxiaua 
of -ICj is at a photon energy less than 0.5 MeV and froa this they 
concluded that the central energy of a grasar line should also be at 
such a photon energy. There are two qualitative effects that cause 
this redshlft. The first one can be seen froa Figure 1. Here, for 
a degenerate pair plasaa of density close to the grasar threshold, 
the naxlmun of occurs very close to 0.5 MeV. But because of 

the steep slope of K^, the aaxlnun of -Kj Is shifted to 0.42 


« 2.2 X lO’® co"3 


«ca*<) 
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MeV. This redshlft, therefore, is a direct coaeequence of the 
effect of the de^eneracf on Cocpton ecsttcring end the high Femi 
eoergr ii^Iied bj the high threshold deoslty* 

Tbs second effect is independent of the degeneracy. It is 
caused by either a high tenperature qj. ^ ^igh dcnalty 

(Feral energy both of uhich broaden and blueohift the 

eoisoivity^* Then froa Eq. (10), division by csoves the peak of 

to an energy <0.5 MeV. In the nondegenerate case, the peck of 
essentially coincides with that of since is not a 

strong function of photon energy. In the degenerate case, with the 
pair density mch higher than the threshold density, the chape of Eq 
does not affect each the position of the tsaxicua of -Ej- 

Varoa^ proposed a graoar nodel based on degenerate pairs %rith 
Feml energy «sc^ for which the peck of would be very close to 

O.Sll t^V. Uen/ever, Coi^ton scattering was ignored in this codel. 

If it were taken into account, then the relatively low pair density 
corresponding to this Fenai energy would Inply that + Kr >0 at 
all energies. Thus grnoar action cannot occur in this oodcl. 

As pointed out in the Introduction, the observed annihilation 
lines in gaaaa-ray bursts are at M).43 KeV. If coce of these lines 
are indeed fron grasar sources, the ieplied gravitational redahift 
would be ouch ecaller than that Icplied by spontaneous, nonacplificd 
annihilation. The potential cidotcnce of graear sources, therefore, 
allows garr!3a*ray bursts to be produced on neutron stare of scalier 
oass and larger radius than previously conjectured, or even objects 
other than neutron stars. 
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